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ABSTRACT 


Surface  mixed  layer  properties  off  Northern  California 
(NOCAL)  were  analyzed  statistically  and  numerically.  The 
observations  were  acquired  on  three  cruises  as  part  of  the 
Pilot  Ocean  Prediction  Study  of  the  California  Current 
eddies  centered  ca.  37  to  39N,  125  to  12 7W  during  March  and 
August  1982.  Mixed  layer  depth,  averaging  33+\4m,  had  a 
horizontal  correlation  scale  of  no  more  than  35km,  which  has 
significance  for  relating  thermal  structure  information  from 
individual  temperature  profiles  to  that  of  Fleet  Numerical 
Oceanography  Center’s  (FNOC)  analyses  based  on  a  grid  length 
of  approximately  300km.  Simulations  and  sensitivity  tests 
were  made  using  the  Garwood  bulk  mixed  layer  model  and  the 
Mellor  Level-2.5  diffusion  model  with  the  initial  and 
boundary  conditions  acquired  at  sea  and  from  FNOC.  Upper 
ocean  thermal  structure  analyses  and  forecasts  were  also 
obtained  from  the  Navy’s  T0PS/T0PS-E0TS  diffusion  model, 
which  has  since  become  operational  at  FNOC.  Comparisons 
of  observations,  analyses,  and  model  solutions  reveal 
consistent  cooling  and  deepening  by  the  former  two  models 
and  excessive  warming  by  the  latter  model.  These 
significant  differences  are  believed  to  be  related  to 
model  resolution,  model  sensitivity,  oceanic  and  atmospheric 

data  quality,  and  spatial  variability. 
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GLOSSARY  OF  NOTATIONS* 


A  Horizontal  eddy  diffusivity  coefficient 

b  Buoyancy 

B  Mean  buoyancy 

b*  Nondimensional  buoyancy 

Sw  Turbulent  buoyancy  flux  (cm**2/s**3) 

cD  Drag  coefficient 

c,Cp  Specific  heat  of  seawater  (cal/gm-C) 

C  Complex  mean  momentum 

cw  Flux  of  complex  momentum 

D  Damping  coefficient  for  inertial  oscillations 

dz  Depth  increment 

e  Total  turbulent  kinetic  energy  (TKE) 

Eq  TKE  flux  from  atmosphere  to  ocean  produced  by 

breaking  waves 

E^  Downward  TKE  flux  at  base  of  homogeneous  layer 
f  Coriolis  parameter  (sec**-l) 

F  Downward  flux  of  solar  radiation 

g  Acceleration  due  to  gravity  (m/s** 2) 

G  Rate  of  energy  production  due  to  mean  shear 

h  Mixed  layer  depth 

k  Wavenumber  of  interfacial  disturbance 

Kh  Vertical  eddy  diffusion  coefficient  for  heat 

H  and  salinity 

Vertical  eddy  diffusion  coefficient  for  momentum 
L,fc  Obukhov  length  scale 
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m, -me  Universal  constants  computed  and  used  in  Garwood  model 

lb 

ML  Mixed  layer 

MLD  Mixed  layer  depth 

p  Pressure 

q  Square  root  of  twice  the  turbulent  kinetic  energy 

Qq  Heat  flux  through  sea  surface 

Heat  flux  through  base  of  ML 
Q(B)  Back  radiation  (cal/sec*cm*ft2) 

Q(C)  Sensible  heat  flux  (cal/sec*cm**2 ) 

Q(E)  Latent  heat  flux  (cal/sec*cm**2 ) 

Q(N)  Net  surface  heat  flux  (cal/sec*cm**2 ) 

Q(S)  Clear  sky  radiation  (cal/sec*cm**2 ) 

Rf  Flux  Richardson  number 

Ri  Gradient  Richardson  number 

Rq  Rossby  number 

S  Salinity 

Sp  Stability  function  for  vertical  eddy  fluxes  of  heat 
and  salinity 

S^j  Stability  function  for  vertical  eddy  fluxes  of  momentum 
SST  Sea  surface  temperature  (C) 
s'w1  Vertical  eddy  flux  of  salinity 
t  Time 

Tq  Temperature  (C)  corresponding  with  sea  surface 

T  Temperature  (C) 

TKE  Turbulent  kinetic  energy 

T *wv  Turbulent  temperature  flux  (mwatts/cm**2) 


u. 


U10,W10 

u,v 

U,V 


u’w’ ,vTw’ 


w 

x,y 


P  » 0 

w  ’  vo 


p 

u 


<> 


o 


Friction  velocity 

Wind  speed  at  10m  (m/sec) 

x  and  y  components  of  current  velocity 

x  and  y  components  of  mean  current  velocity 

Momentum  fluxes 

z  component  of  current  velocity 

Grid-referenced  horizontal  corrdinates , 
positive  toward  east  and  north,  respectively 

Vertical  coordinate,  positive  upward  from  sea 
surface 

Reference  density  for  water 

Density  of  seawater 

Molecular  diffusivity  coefficient 

Expansion  coefficient  for  heat 

Expansion  coefficient  for  salt 

Vertical  average 

Ensemble  mean 

Departure  from  mean 

Surface  wind  stress  ( dvnes/cm**2 ) 

Dissipation  time  scale 
Time  scale  of  TKE  transport 
Dissipation 


*  Attempts  have  been  made  to  maintain  original  notations. 
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I.  INTRODUCTION 


A.  BACKGROUND  AND  OBJECTIVES 

Over  the  past  few  years,  the  ocean  research  community 
and  the  Navy  have  become  increasingly  aware  of  the  need  and 
measures  for  ocean  thermal  structure  modeling.  Not  only  has 
the  ocean  response  to  atmospheric  forcing  become  of  impor¬ 
tance  because  of  its  feedback  to  the  atmosphere,  but  also 
because  of  the  importance  of  ocean  thermal  structure  in 
acoustic  propagation.  For  these  and  other  reasons,  a  grow¬ 
ing  study  has  centered  on  the  upper  ocean-  "that  layer 
bounded  by  the  surface  which  exhibits  strong  seasonal  varia¬ 
tions... most  notable  in  the  local  temperature  profile  of  the 
water  column.”  [Grabowski,  et.al.,  1982]  The  nearsurface 
structure  called  the  mixed-  or  mixing  layer  has  received 
significant  emphasis.  The  Glossary  of  Meteorology  defines 
this  structure  "as  a  surface  layer  of  virtually  isothermal 
water  due  to  wind  induced  turbulent  motion  and/or  free 
convection  which  frequently  exists  above  the  thermocline” . 

To  understand  the  form  and  dynamics  of  the  mixed  layer  and 
the  upper  ocean  thermocline  in  general,  numerous  models  have 
come  into  use. 

The  objective  of  this  research  is  first  of  all  to  obtain 
a  description  of  the  characteristics  of  the  mixed  layer  in  a 
'test  block'  of  the  ocean  (ca.  37  to  39N;  125  to  127W) 


surveyed  during  the  first  and  second  OPTOMA  cruises  (Figs.  1,  2) 
occurring  during  the  second  week  of  March  and  the  first  two 
weeks  in  August  1982,  respectively.  (OPTOMA-  Ocean  Prediction 
Through  Observations,  Modeling,  and  Analysis,  -is  a  joint 
Harvard/NPS  project  ’intended  to  acquire  field  data  to 
characterize  synoptic  scale  eddies  over  a  domain  in  the 
California  Current  off  Northern  California,  and  to  'set  up'  an 
eddy-resolving,  statistical/dynamical,  limited-domain,  open 
boundary  numerical  ocean  prediction  model'.)  The  next 
objective  was  to  evaluate  the  thermal  structure  analysis  of 
the  Navy's  TOPS  model  over  the  second  cruise  period  in  com¬ 
parison  with  the  Mellon  level-2.5  model  and  the  Garwood 
(1976,  1977)  bulk  mixed  layer  model.  Finally,  the  results  of 
the  three  basically  single  point  models  were  compared  with  the 
observations.  (Note:  TOPS  is  often  considered  a  quasi- 
three-dimensional  model,  because  it  includes  advection  by 
climatological  geos trophic  currents  and  wind  drift.) 

The  idea  that  the  vertical  distribution  of  temperature, 
salinity,  and  current  in  the  upper  ocean  layers  generally 
appears  to  be  governed  largely  by  the  vertical  heat,  momentum, 
and  salt  fluxes  imposed  by  local  air-sea  transfer  [cf. 

Camp  and  Elsberry,  1978]  has  led  to  one-dimensional  time 
dependent  models  which  assume  horizontal  homogeniety. 

[Klein,  1980]  The  reliance  upon  this  assumption  varies 
among  models  as  well  as  regions  of  applicability. 
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B.  LITERATURE  REVIEW 


Techniques  for  modeling  vary  significantly  in  the  degree 
to  which  they  attempt  to  represent  the  upper  ocean  features 
and  dynamics.  The  existing  models  may  be  categorized  as 
belonging  to  one  of  two  broad  classes: 

a.  Bulk  Mixed-Layer  Models 

b.  Diffusion  Models 

Each  class  of  models  will  be  described  with  some  reference 
to  its  components.  No  attempt,  however,  will  be  made  to  cover 
all  of  the  models  in  great  detail  or  to  describe  each  of  the 
existing  models  or  model  modifications.  The  only  exceptions 
to  this  will  be  the  Garwood  (1977)  and  TOPS  [Mellor  and 
Durbin,  1975]  models  which  are  of  direct  interest  in  this 
study.  The  structure  of  the  Mellor  Level-2.5  model  will  not 
be  specifically  discussed,  because  it  is  very  similar  to  the 
TOPS  (Mellor  Level-2)  model. 

1.  Bulk  Models 

The  first  category  of  models,  the  bulk  models,  are 
often  referred  to  as  slab  models.  The  reason  for  this  is 
that  Pollard  and  Millard  (1970)  and  Halpern  (1974)  showed, 
with  the  occurrence  of  strong  wind-forced  vertical  mixing, 
the  mixed  layer  responds  at  the  inertial  frequency  as  a 
rigid  slab.  This  description  may  be  misleading,  for  it 
implies  that  T,  S,  and  V  are  uniform  throughout  the  layer. 
(Glossary  of  Tents)  Such,  however,  is  not  the  case.  Instead, 
the  models  are  based  on  the  assumption  that  values  of  the 


aforementioned  parameters  are  ’quasi-uniform’  within  a 
'well-mixed  layer’.  [Klein,  1980]  Provided  the  assumption 
holds  (as  it  may  away  from  fronts  and  eddies),  the  values 
can  be  'lumped'  into  integrals.  [Niiler  and  Kraus,  1977; 
Zilitinkevitch  et.al,  1979]  That  is,  the  'bulk'  models 
predict  integral  (or  average  mixed  layer)  values  for  all 
variables.  (Garwood,  personal  communication) -hence ,  the 
perhaps  'better'  term  -'integral  model'  is  frequently  used. 

This  type  of  modeling  requires  an  a  priori  assumption 
that  a  mixed  layer  exists.  [Marchuk,  et.al.,  1977]  To 
describe  the  layer  requires  only  knowledge  of  the  fluxes 
at  the  surface  and  the  base  of  the  mixing  layer  (i.e., 
surface  and  entrainment  fluxes).  The  mathematics  of  the 
modeling  is,  therefore,  greatly  simplified  from  partial 
differential  equations  to  ordinary  differential  equations. 

Basically,  mixed  layer  modeling  is  founded  on  two 
fundamental  hypotheses.  The  first  is  that  the  mixed  layer  is 
formed  as  a  result  of  the  upper  ocean  response  to  forced  and 
free  convection.  Atmospheric  forcing  in  the  form  of  applied 
wind  stress  converts  mean  kinetic  and  potential  energy  to 
turbulent  kinetic  energy  within  the  upper  ocean.  This 
turbulent  motion,  in  turn,  entrains  below -layer  water  into  the 
layer,  thereby  changing  the  layer's  characteristics.  Strong 
solar  heating  works  against  the  forced  convection  to 
increase  the  potential  energy,  and  form  a  shallow  turbulent 
boundary  layer.  This  heating  in  effect  alters  only  the  near 
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surface  water,  for  invoking  the  Second  Law  of  Thermodynamics, 
water  cannot  be  unmixed.  Thus,  solar  heating  results  in  a 
new  mixed  layer  in  addition  to  the  'old  mixed  layer'. 

(Fig.  3)  This  distinction  between  new  and  old  mixed  layers 
will  be  discussed  and  utilized  in  the  following  sections. 

The  second  hypothesis  of  importance,  a  topic  of  much 
controversy,  invokes  the  use  of  some  form  of  the  mechanical 
energy  budget  as  essential  in  the  closure  scheme  of  the  more 
recent  models.  The  disagreement,  and  consequently  a  differ¬ 
ence  in  some  models,  lies  in  the  question  as  to  whether  the 
kinetic  energy  of  the  mean  flow  or  the  turbulent  kinetic 
energy  plays  the  main  role  in  driving  the  downward  buoyancy 
flux  at  the  base  of  the  turbulent  boundary  layer,  i.e.,  in 
controlling  the  mixed  layer  dynamics  and  entrainment. 

Thompson  (1976),  Pollard  (1977),  and  Price  (1977)  advocate 
that  the  mean  kinetic  energy  is  converted  directly  to  poten¬ 
tial  energy  by  an  instability  in  the  mean  flow  that  causes 
the  bulk  Richardson  number  to  be  less  than  some  constant  of 
order  one.  Advocates  of  the  turbulent  kinetic  energy 
budget,  following  the  original  argument  of  Kraus  and  Turner 
(1967),  disagree  that  this  will  be  the  dominant  mechanism, 
reasoning  that  turbulence-generated  instabilities  normally 
have  the  capacity  to  erode  the  interface  to  the  extent  that 
they  preclude  the  existence  of  mean  flow  instabilities. 

C Garwood,  1979] 
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Figure  3.  Temperature  versus  depth  profile  with  old 
and  new  MLD's;  h  =new  mixed  layer  depth, 
h„=old  mixed  layer  depth  [Husby  and 
Nelson,  1982] 


Differences  in  the  views  by  no  means  end  there;  for 
example,  there  is  no  unanimous  approach  to  parameterization 
of  the  kinetic  energy  budget.  The  gross  parameterization 
approach  is  commonly  used  in  closure  schemes.  However, 
second  order  closure  schemes  have  been  developed,  for  example 
Launder  (19  7  5).  Such  closure  schemes,  though  applied  in  only 
one  of  the  bulk  models  -  Garwood  (1976,  1977)-are  receiving 
increasing  use.  Garwood's  scheme  will  be  described  in  a 
later  section. 

2 .  Summary  of  Bulk  Models 

The  structure  of  all  simple  one-dimensional  bulk 
models  may  be  described  in  terms  of  two  fundamental  equations 
the  vertically  integrated  heat  equation  and  the  vertically 
integrated  turbulent  kinetic  energy  (TKE)  budget. 

(1)  h(3T0/3t)  =  Qq  -  0h 

A  B  C 

h  h  , 

(2)  P  =  /  Gdz+EQ-Eh+S(o/  Fdz-  |hF|z=h) 

D  E  F  G  H1  H2 

h  h  ae 

-  /  edz  -  /  II-  dz. 
o  o  3T 

I  J 
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where  indicated  terms  are: 

A:  Mixed  layer  depth  (MLD)  times  the  local  time  rate  of 

change  of  mixed  layer  temperature  (MLT). 

B:  Net  heat  flux  through  the  sea  surface. 

C:  Net  heat  flux  through  the  base  of  the  ML. 

D:  P=  ieh(Q  +  Q,  ) 

Vertically  integrated  generation  or  dissipation  of 
TKE  by  buoyancy  forces. 

E:  Rate  of  TKE  production  due  to  transfer  from  mean 

current  shear. 

F:  E  is  the  downward  TKE  flux  at  the  top  of  the  ML 

o 

(i.e.,  from  the  atmosphere  to  the  ocean  and 
produced  by  wave  breaking  at  the  sea  surface). 

G:  E^  is  the  downward  turbulent  energy  flux  at  the  base 

of  the  ML. 

and  :  Effect  of  penetrative  solar  radiation. 

I:  Net  turbulent  layer  dissipation  over  the  ML. 

J:  Net  local  time  rate  of  change  of  TKE  over  the  ML. 

The  main  differences  between  the  various  bulk  models 
are  in  the  definitions  of  the  heat  flux,  Q,  and  the  turbu¬ 
lent  entrainment  across  the  interface.  Some  of  the  main 
features  of  several  bulk  models  are  listed  below. 

Kraus  and  Turner  (KT)  (1967);  KT  developed  the  first  bulk 
model  capable  of  addressing  the  time  evolution  of  a  mixed 
layer/seasonal  thermocline  system;  they  were  the  first  to 
propose  a  TKE  budget  to  close  the  system.  Kraus  and  Turner 


assumed  a  steady  state  turbulent  regime  and  neglected 
dissipation.  They  parameterized  penetrative  solar  radiation 
with  an  exponential  decay  law  and  shear  production  in  terms 
of  the  frictional  velocity. 

Denman  (DN)  (1973):  DN  rederived  the  KT  model ,  following 
Kraus'  advice  that  the  rate  of  increase  of  potential  energy 
by  mixing  should  be  a  small  constant  fraction,  m,  of  the 
rate  at  which  turbulent  kinetic  energy  is  transferred  down¬ 
ward  by  the  wind  stress.  (He  used  m=.0012)  Denman  also 
included  absorption  of  solar  radiation  below  the  mixed 
layer.  He  showed  that  for  sufficiently  long  model  runs,  the 
eventual  depth  of  the  mixed  layer  depended  as  much  upon  the 
temperature  gradient  below  the  mixed  layer  as  upon  the 
generation  of  turbulent  kinetic  energy.  Unlike  KT ,  DN 
included  some  turbulent  dissipation. 

Pollard,  et.al.  (PRT)  (1973):  As  DN  and  KT,  PRT  assumed  a 
steady  state  regime.  The  rate  of  work  done  by  the  wind 
stress  on  the  mean  motion  was  equated  with  the  rate  of 
increase  of  the  energy  of  the  mean  field.  The  mixed  layer 
deepened  rapidly  to  a  maximum  within  half  of  an  inertial 
period  from  the  time  of  the  instantaneous  increase  in  the 
wind  speed  (which  then  was  held  constant).  Forty  percent 
of  the  incoming  energy  was  found  to  directly  increase  the 
mean  kinetic  energy  while  the  rate  of  change  of  the  mean 
potential  energy  was  three  times  smaller  than  that  of  the 
mean  kinetic  energy.  This  result  reinforced  their 


28 


hypothesis  that  the  shear  instability  of  non-steady  ’wind 
drift  currents'  is  the  main  mechanism  in  turbulence 
production.  They  neglected  turbulent  dissipation  and 
penetrative  solar  radiation. 

Niiler  (NI  ’75)  (1975):  MI  recognized  that  instantaneous 
sharp  wind  changes  were  not  realistic,  and  that  in  a 
steady  regime,  the  surface  energy  fluxes  from  the  breaking 
waves  could  be  an  order  of  magnitude  greater  than  the  energy 
fluxes  from  the  wind  to  the  currents.  Thus,  NI  extended  the 
PRT  model  so  that  at  the  beginning  and  final  stages  of  a 
model  run,  TKE  generation  through  surface  fluxes  was  most 
important,  while  at  intermediate  times  the  drift  current 
velocity  shear  was  dominant.  As  with  the  PRT,  NI'75 
assumed  a  steady  state,  no  penetrative  solar  radiation,  and 
no  energy  flux  through  the  base  of  the  mixed  layer. 

Turbulent  dissipation  was  parameterized  as  a  function  of 
the  frictional  velocity. 

Kim  (KIM)  (1976)  :  KIM  proposed  a  model  more  complex  than  KT 
which  included  TKE  storage  as  well  as  the  same  parameteri¬ 
zation  for  penetrative  solar  radiation  as  in  DN.  An 
additional  feature  of  importance  was  the  use  of  a  background, 
depth-dependent  turbulent  dissipation  rate  to  avoid  infinite 
deepening. 

DeSzoeke  and  Rhines  (DR)  (1976);  They  used  the  NI  model 
with  the  addition  of  a  TKE  storage  term  (term  J)  to  show 
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that  the  KT  and  PRT  models  yielded  results  that  are 
assymptotes  for  a  single  model  equation. 

Gill  and  Turner  (GT)  (1976):  They  simulated  oceanic 
seasonal  cycles,  and  showed  that  the  cyclical  behavior  of 
the  potential  energy  content,  heat  content,  and  the 
surface  temperature  could  be  recovered  if  all  or  some 
portion  of  the  potential  energy  released  during  the  runs 
was  dissipated.  They  neglected  terms  G,  H^,  H2,  and  J. 

Niiler  (NI  *77)  (1977);  This  model  is  similar  to  NI  '75. 
Additions  include  a  TKE  storage  term  (term  J)  and  background 
dissipation  identical  to  that  employed  by  Kim. 

Niiler  and  Kraus  (NK)  (1977):  The  NK  model  is  a  conglomerate 
of  several  of  the  models  described  above.  For  example,  the 
radiation  term  is  the  same  as  that  used  in  the  KT  model; 
terms  E  and  F  are  similar  to  those  used  previously  by 
Niiler  with  the  only  difference  being  in  the  coef f icients ; 
and  the  turbulent  dissipation  term  is  similar  to  that  used 
by  Resnyansky  (1975)  and  the  shear  production  term  of  the  GT 
model.  The  NK  model  does  not  include  the  storage  term,  J. 
Unlike  previous  models,  it  includes  the  downward  flux  of  TKE 
through  the  mixed  layer  base  (term  G) . 

Garwood  (GWD)  (1977):  (to  be  discussed  in  Section  3) 


TABLE  I** 


SUMMARY  OF  TERMS  USED  IN  BULK  MODELS 
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*  Appears  explicitly 

+  Appears  implicitly,  i.e.,  in  one  of  the  other  terms 
**  Modified  (Zilitinkevich,  et . al .  ,  1979) 


3.  Garwood  Model 


The  Garwood  (1977)  bulk  mixed  layer  model  is 
composed  of  a  closed  system  of  seven  equations: 

The  entrainment  flux  equation: 


(3) 


~7  j  - 
mu<w  >  ‘  <E> 

-  Ew(-h)  =  - c - 


The  horizontal  component  of  the  turbulent  kinetic  energy 
equation: 


(4) 
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The  vertical  component  of  the  turbulent  kinetic  energy 
equation: 


(5)  ^i_(h<w^>)=  j  hbw(-h)-  ^  hu^b^+m^  (  <E>-3  ( ) )  <E> 


m.  ,/2  nir 
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fh) <E> . 


The  mean  buoyancy  and  complex  mean  momentum  equations: 


(6) 
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(7)  h  - -  =  cw(-h)-cw(0)-if <C>h. 

3  c 
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The  jump  conditions  (relating  entrainment  fluxes  to  mean 
buoyancy,  rate  of  deepening,  and  mean  momentum): 

(8)  -5w( -h )  =  AB  f—p  . 

0  t 

(9)  -cw(-h)=AC  , 

Depending  on  which  version  of  the  model  is  used, 
one  needs  either  to  provide  surface  momentum  and  buoyancy 
fluxes,  mean  buoyancy  and  momentum  below  the  mixed  layer 
along  with  the  radiation  absorption  or  else  the  observed 
variables  from  which  these  fluxes  can  be  computed.  Closure 
is  achieved  with  the  bulk  buoyancy  and  momentum  equations 
and  the  vertically-integrated  individual  turbulent  kinetic 
energy  components. 

A  number  of  significant  differences  exist  between 
this  model  and  other  bulk  mixed  layer  models  that  may  not  be 
obvious  from  the  preceding  summarization  or  the  matrix  of 
terms  describing  which  terms  are  included  in  each  model. 

Some  of  these  will  be  made  explicit  in  the  following 
discussion  of  the  mechanics  of,  and  ideas  beind,  the  Garwood 
mode  1 . 

According  to  Garwood  (1977),  the  basis  for  the 
initial  interface  destabilization  and  the  resulting  entrainment 


is  the  shear  across  the  interface  (base  of  the  mixed  layer). 
This  shear  is  hypothesized  by  Garwood  to  be  the  result  of 
local  turbulent  eddies  (wind  generated  turbulent  kinetic 
energy  available  for  mixing  is  explicitly  parameterized  as 
dependent  on  h/L)  with  the  mean  flow  only  making  a  minor 
contribution  to  the  achievement  of  this  critical  shear 
value.  (The  instability  related  to  this  criterion  is  often 
called  a  Kelvin-Helmholtz  instability  or  a  Benjamin  (1963) 
Class  C  instability. )  By  using  the  familiar  equation  for 
the  flux  Richardson  number: 

(10)  Rf =Bw/ ( uwSU/3 z+vw3 V/3 z )  ,  and 

taking  into  account  the  jump  conditions  (8,  9),  shear 
production  can  be  shown  to  be  only  a  fixed  fraction  of  the 
buoyant  damping  in  the  entrainment  zone.  This  zone  may  in 
fact  have  a  Rf>  1,  and  still  possess  enough  turbulent 
kinetic  energy  for  entrainment  to  proceed.  Garwood's 
reasoning,  in  contrast  to  that  of  Pollard, et .al .  (1973), 
is  that  the  mean  shear  production  is  only  a  secondary 
energy  source  for  mixing.  This  secondary  source  of  energy 
can  only  be  made  available  by  entrainment  which  is 
initiated  by  another  source.  This  other  and  more 
significant  source  is  the  convergence  of  the  flux  of 
turbulent  energy  from  above  by  the  term: 

(11)  C w ( ■«-  ♦  -2 — )] 

3Z  i  0q 
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Therefore,  the  critical  number  in  the  Garwood  model  is  not 
the  flux  Richardson  number,  but  instead  the  ratio  of  the 
buoyancy  flux  to  the  convergence  of  turbulent  energy,  P: 


P=Ew/|— [w(£-  +  -£-)] 
3  z  2  o  _ 


The  t ime ,  t  , 
e 


.  ~~7  -1/2 
t  =  a,h<w  > 
e  1 


required  to  transport  the  available  turbulent  energy,  <E> , 
to  the  entrainment  zone  is  important  in  relation  to  dissipa¬ 
tion,  and  to  the  idea  that  not  all  available  TKE  contributes 
to  increasing  the  potential  energy  of  the  system.  According 
to  Tennekes  and  Lumley  (1972),  the  viscous  dissipation  of 
geophysical  flows  with  large  Reynolds  numbers  is  proportional 
to  the  reciprocal  of  the  time  scale  of  the  largest  eddies. 
Garwood  developed  the  idea  that  for  the  deeper  mixed  layers 
(Ro~l),  the  planetary  rotation  introduces  a  second  time  scale 
of  importance:  (1/f).  Not  only  did  the  inclusion  of  this 
second  time  scale  enable  Garwood  to  more  explicitly 
parameterize  viscous  dissipation  in  relation  to  the  local 
Rossby  number,  but  the  very  role  of  the  time  scale  in  the 
entrainment  equation  enabled  Garwood’s  model  to  simulate  a 
wider  range  of  conditions,  i.e.,  diurnal  and  annual  ranges 
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of  mixed  layer  stability.  (The  PRT  model  can  simulate  a  wide 
range  of  changes,  though  the  deepening  of  the  mixed  layer  is 
in  accordance  only  with  the  mean  shear  production.) 

Another  difference  over  most  previous  models  was  in 
Garwood's  method  of  using  the  turbulent  kinetic  energy 
equation.  The  model  uses  the  component  equations  (4,  5) 
vice  the  typical  use  of  the  total  equation.  This  facilitates 
analyzing  the  retreat  or  shallowing  of  the  turbulent  boundary 
layer  which  depends  on  only  the  vertical  component  of  turbu¬ 
lence  being  inadequate  to  transport  heat,  momentum,  and  TKE 
to  the  earlier  depth  of  mixing.  (Fig.  3)  Thus,  using  the 
equation  in  component  form  allows  for  a  more  explicit  and 
satisfactory  treatment  of  the  mixing  process. 

4.  Diffusion  Models 

The  second  category  of  mixed  layer  models,  the  dif¬ 
fusion  models,  has  received  less  attention  from  the  ocean¬ 
ographic  community  than  the  simpler  bulk  models.  Yet,  the 
interest  that  exists  has  grown  greatly  with  the  increasingly 
powerful  computers  and  the  development  of  methods  for  solving 
nonlinear  partial  differential  equations.  [Marchuk,  et.al 
1977]  The  diffusion  models  are  both  complex  to  interpret 
and  to  program.  They  require  the  numerical  solution  of 
partial  differential  equations  of  temperature,  salinity,  and 
momentum  in  time  and  depth  for  closely  spaced  vertical 
arrays  of  gridpoints.  Unlike  bulk  models,  the  'gridpoint 
models'  do  not  require  the  a  priori  assumption  of  the 
existence  of  a  mixed  layer.  Mixed  layer  depths,  are,  on  the 
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contrary,  diagnosed  from  the  examination  of  the  turbulence 
flux  predictions  as  well  as  the  predicted  temperatare  and 
salinity  profiles. 

Mathematically,  the  ’mixed  layer'  in  a  diffusion 
model’s  vertical  profile  is  treated  in  the  same  manner  as 
is  the  thermocline.  This  allows  for  a  more  detailed  analysis 
of  the  vertical  structure  (mixed  layers,  thermoclines ,  etc.), 
the  mixing  processes  throughout  the  water  column,  and  a  more 
complete  simulation  of  the  vertical  fluxes  of  temperature, 
salinity,  and  momentum.  The  latter  is  a  basic  concern  in 
the  diffusion  models  which  attempt  to  represent  the  turbulent 
fluxes  w'u',  w'T ' ,  and  w' S v  by  algebraic  or  differential 
equations  in  terms  of  the  gradients  of  mean  temperature, 
salinity,  and  momentum.  These  fluxes  in  conjunction  with 
the  necessary  boundary  conditions  for  T,  S  and  u  provide  a 
solution  for  the  basic  conservation  equations  of  temperature, 
salinity,  and  momentum.  Thus,  the  system  of  turbulence  and 
thermodynamic  equations  are  closed  in  order  to  solve  for 
the  vertical  distributions  of  thermodynamic  and  turbulence 
variables.  [Grabowski,  1980;  Kondo,  et.al.,  1978] 

Diffusion  models  may  be  further  categorized  by 
complexity,  i.e.,  a)  no-equation  models;  b)  one-equation  or 
k-models;  c)  two  equation  or  k-l-models;  and  d)  second 
order  closure  models. 

A  no-equation  model  is  based  on  the  turbulence 
transport  being  parameterized  by  a  buoyancy  influenced  eddy 
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viscosity  coefficient.  The  one-equation  or  k-models  are 
similar  to  the  no-equation  models  though  they  include  a 
relationship  between  eddy  diffusion  and  the  local  turbulence. 
The  two-equation  or  k-l-models  include  a  transport  differen¬ 
tial  equation.  The  second  order  closure  models,  first 
proposed  by  Chou  [1945],  are  based  on  transport  equations 
for  all  of  the  Reynolds  tresses  derived  from  the  Navier 
Stokes  equations.  They  also  include  turbulent  heat  and 
salinity  fluxes  and  the  variance  derived  from  the  latter  in 
conjunction  with  the  heat  conservation  equation. 

A  detailed  profile  is  sought  to  represent  the 
oceanic  vertical  structure  as  completely  as  possible.  For 
example,  Halpern  (1976)  among  others,  noted  the  existence  of 
shear  in  the  top  of  the  water  column,  a  characteristic 
assumed  net  to  exist  in  bulk  models.  Subsequently,  the 
necessity  of  using  numerical  models  to  determine  a  more 
realistic  distribution  of  variables  such  as  the  velocity, 
temperature,  turbulence  intensity,  dissipation,  and  the  like 
was  recognized  [Phillips,  1977], 

The  first  ocean  diffusion  model  was  implemented  by 
Ekman  in  1905.  It  was  a  type  A  model,  that  is,  it  was  a 
no-equation  model  that  neglected  buoyancy  and  salinity.  A 
notable  result  from  this  model  was  the  famous  Ekman  spiral. 
Because  of  its  neglect  of  density  stratification,  this  model 
will  not  be  discussed  further.  The  major  difference  between 
the  diffusion  models,  beyond  the  number  of  equations 
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employed,  lies  in  whether  they  are  based  on  observations, 
dimensional  analysis,  and  simple  stability  arguments  or  on 
some  form  of  the  TKE  energy  equation.  Of  the  five  models  to 
be  briefly  discussed  below,  the  first  two  are  of  the  former 
type,  while  the  other  three  are  of  the  latter  type. 

5 .  Summary  of  Diffusion  Models 

Munk  and  Anderson  (MA):  1 ? 4 9 :  The  MA  model  is  a 

no-equation  model  wherein  the  eddy  diffusion  coefficients 
for  heat  and  momentum  are  functions  of  the  gradient 
Richardson  number.  Mo  salinity  is  included  in  this  model. 

Kondo,  Sasano,  and  Ishii  (KSI):  1979:  The  main  fea¬ 
ture  of  the  SKI  model  is  that  constant  flux  formulations  are 
applied  to  regions  of  varying  fluxes.  Formulation  of 
diffusion  coefficients  based  on  studies  of  the  atmospheric 
boundary  layer  are  cast  in  Monin  Obukhov  formalism  and 
aDDlied  to  the  upper  oceans.  Eddy  diffusion  coefficients 
for  salt  are  assumed  to  be  the  same  as  those  for  heat. 

Vager  and  Zilitinkevich  (VZ):  1958:  This  model  is 

a  one-equation  or  k-model  based  on  Prandtl's  extension  of 
the  mixing  length  theory  for  homogeneous  flow.  Salinity  is 
ignored . 

Mellor  and  Yamada  (MY):  1974;  Mellor  and  Yamada 

formulated  a  sequence  of  models  with  varying  degrees  of 
complexity  using  second  order  closure.  This  discussion  will 
deal  with  the  level-2  model,  because  the  TOPS/TOPS-EOTS 
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model  (to  be  discussed  in  the  next  section)  is  based  on  the 


level-2  MY  turbulence  closure  scheme,  and  because  it  is  simpler 
than  the  level-2.5  model.  (The  level-2.5  MY  model  was  also 
used  though  its  mechanics  are  quite  similar  to  those  of  the 
level-2  MY  model.  The  main  difference  is  the  inclusion  of 
additional  complexity  as,  for  instance,  the  retention  of  the 
TKE  storage  term  in  the  level-2.5  MY  model.)  (To  avoid 
repetition,  TOPS/TOPS-EOTS  will  not  be  mentioned  in  this 
section  except  to  say  that  the  ’Clancy  version'  employed 
includes  salinity,  solar  radiation,  and  horizontal  advection 
by  instantaneous  wind  drift  and  climatologically  averaged 
geostrophic  currents.)  The  level-2  MY  model  neglects 
horizontal  advection  and  diffusion  terms.  The  TKE  budget  is 
consequently  a  balance  between  the  production  and  dissipation 
of  TKE  and  the  conversion  to  potential  energy  as  seen  in 
Table  II.  A.  similar  balance  exists  between  the  turbulent 
heat  flux,  the  temperature  variance,  and  the  individual 
stress  components.  Unlike  the  other  two  models  described  in 
Table  II,  the  MY  level-2  model  neglects  the  turbulent  transport 
term. 

Kochergin,  Klimok,  and  Sukhoruks  (KKS):  1976:  The 
KKS  model  is  a  two  equation  or  a  k-l-model.  Its  TKE  budget 
parameterizes  the  time  rates  of  change  of  TKE  and  the 
turbulence  dissipation  rate. 
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table  II 


SUMMARY  OF  TERMS  USED  IN  DIFFUSION  MODELS 


(Modified: 

Grabowski 

,  et . 

al.  1980) 

TKE  Prod/  Shear 
dissip  Produc- 

rate  tion 

Buoyancy 
Conver¬ 
sion  of 

PE  to  TKE 

Vertical 
Transport  of 
TKE  Work  Done 
against  Pres¬ 
sure  Gradient 

Dissipation 

vz 

* 

* 

A 

* 

0 

MY 

0 

A 

A 

0 

0 

A 

KKS 

A 

A 

A 

* 

0 

A 

includes 

term 

0:  does  not  include  term  explicitly 


6.  TOPS/TOPS-EOTS 

T0PS/T0PS-E0TS  is  the  Navy's  new  real-time  ocean 
thermal  structure  analysis/forecast  system  in  use  at  the 
Fleet  Numerical  Oceanography  Center  (FNOC).  (Operational  as 
of  March  1983)  It  was  developed  by  the  Naval  Ocean 
Research  and  Development  Activity  (NORDA)  as  part  of  the  Navy's 
automated  environmental  prediction  system  (AEPS).  The  goal 
behind  TEOTS  is  to  produce  reliable,  accurate,  real-time 
representations  of  the  ocean  thermal  structure.  The  motivation 
for  the  way  the  model  is  organized  (to  be  discussed  below)  lies 
in  the  available  sources  of  data  and  the  operational  nature 


of  the  model.  Given  that  the  upper  ocean  is  primarily 
atmospherically  forced,  and  the  assessment  that  the  mixed 
layer  depths  have  proved  ’highly  predictable  with  a  variety 
of  models’,  such  as,  the  Denman  (1973)  or  Mellor  and  Durbin 
(MD)  (1975)  models,  NORDA  has  developed  an  operational  model 
based  largely  on  the  use  of  the  rather  well-defined  atmos¬ 
pheric  variables  to  drive  a  model  to  describe  the  relatively 
data  sparse  oceans.  Surface  wind,  solar  heat  fluxes,  and 
precipitation  fields  are  among  those  necessary  to  drive 
""OPS/TEOTS. 

For  vertical  structure,  TOPS  has  seventeen  levels 
from  the  surface  to  500  meters,  with  the  majority  of  the 
levels  placed  within  the  first  100m  to  enhance  the  vertical 
resolution  where  most  needed.  (Appendix  B)  The  horizontal 
grid  is  typical  of  many  operational  atmospheric  models, 
i.e.,  381  km  at  60N.  The  input  fields  on  this  grid  are 
smoothed  spatially  by  the  TOPS  modified  EOTS.  Spatial 
averaging  of  the  fields  also  smooths  the  temporal  varability 
of  the  fields.  Of  course,  horizontal  features  smaller  than 
the  spacing  cannot  be  resolved  so  that  the  capabilities  are 
limited  in  enclosed  seas,  in  the  vicinity  of  fronts  and 
eddies,  and  near  boundary  currents.  [Clancy  and  Poliak,  1982] 
This  last  point  is  important  to  keep  in  mind  for  the 
simulation  comparisons  in  Chapter  IV. 
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As  its  name  implies,  TOPS/TOPS-EOTS  is  made  up  of 
two  components:  TOPS  and  TOPS-EOTS.  They  are  coupled  in  a 
cyclical  fashion. 

TOPS-EOTS  is  an  improved,  updated  version  of  the 
Navy's  conventional  objective  analysis  component,  EOTS 
(Expanded  Ocean  Thermal  Structure).  It  supplies  TOPS  with 
updated  (real-time)  initial  conditions  for  a  twenty-four 
hour  forecast.  Its  improvement  over  conventional  EOTS  lies 
in  the  fact  that  the  TOPS  forecast  is  fed  back  into  TOPS-EOTS 
as  a  first  guess  field  for  the  following  day's  analysis.  This 
not  only  reduces  the  noise  resulting  from  the  introduction  of 
data  into  the  analysis,  and  makes  the  analysis  dynamically 
and  thermodynamically  consistent  with  the  atmospheric 
forcing,  but  it  also  allows  the  model  to  be  run  for  an 
extended  period  of  time  within  data-sparse  regions  prior  to 
being  updated,  provided  atmospheric  forcing  data  are  available . 

TOPS  (Thermodynamic  Ocean  Prediction  System),  the 
forecast  component,  is  a  thermodynamic  model  for  vertical 
structure  which  is  forced  by  surface  fluxes  supplied  by 
NOGAPS  (Naval  Operational  Global  Atmospheric  Prediction 
System).  It  uses  the  MY  level-2  turbulence  scheme  previously 
described  with  the  added  effects  of  advection  by 
instantaneous  wind  drift  and  climatologically  averaged 
geostrophic  currents. 


43 


The  prognostic  equations  used  by  the  model  are  the 


conservation  equations  of  temperature,  salinity  and  momentum 


(14) 


aT  _  i_ 

at  ~  az 


(=wTTr  + 


+  JL  II 

pwc  az 


a_ 

ax 


(u  T)  -  f-(v  7)  -  f-(w  T)  +  A( 
a  3y  a  az  a 


+ 


3y 


(15)  |f 
3  T 


a  z 


( -w'  S 1  + 


as,. 

',32) 


a 

ax 


(u  S)  -  |-(v  5)  -  f-(w  5 
a  3y  a  3  z  a 


)  +  A( 


O 


(16) 


au  _ 
at 


fv  +  (-w’uT 

3  Z 


+ 


u 

az 


-Du 


a  v  _ 
at 


-fu  +  ~  ( -W ’ V '  + 
az 


-Dv 


where : 

'D'  is  a  damping  coefficient  representing  the  drag  force  by 
radiational  stress  at  the  base  of  the  mixed  layer  associated 
with  the  propagation  of  internal  waves  downward  and  away 
from  the  wind-forced  region.  ’A’  is  the  background 


horizontal  eddy  diffusion  associated  with  the  intermittent 
breaking  of  internal  waves.  [Clancy,  et.al.,  1982] 

There  are  no  horizontal  pressure  gradients  in  the 
last  two  equations  so  that  the  u  and  v  represent  wind  drift 
current.  The  geostrophic  advection  is  provided  to  the  u 
and  v^  terms  in  the  T  and  S  equations. 

Within  the  MY  level-2  closure  scheme,  the  vertical 
eddy  fluxes  of  T,  S,  and  momentum  are  parameterized 
(Mellor  and  Durbin  (1975);  Clancy  and  Poliak  (1983))  as 
follows : 
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where  K„  and  K„  are  eddy  diffusion  coefficients;  i  is  a 
M  H 

turbulence  length  scale;  q=sqrt ( ( 2 )*TKE) ;  and  are 
functions  of  the  bulk  Richardson  number: 
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(20)  Ri  s 


There  is  an  implied  cutoff  Richardson  number  such 
that  one  can  determine  the  three  empirical  constants.  The 
TKE  budget  used  to  obtain  q  is: 


(21)  tqSH[(if)2  *  (||)2]  ♦  tqSH<8-  §£)  -  fa  -  0 


The  solar  radiation  flux  is  parameterized  by  using  an 
extinction  profile  for  the  most  common  water  type.  (For 
the  level-2.5  model,  the  solar  radiation  flux  is  applied 
at  the  surface.) 

As  mentioned  previously,  boundary  conditions  are 
mainly  provided  by  the  FNOC  atmospheric  models.  XBT  profiles 
are  included  in  the  analysis  when  they  are  available.  (The 
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ocean  data  used  in  the  TOPS  runs  will  be  briefly  mentioned 
in  later  sections  when  the  model  results  are  discussed.) 


7 .  Verification  Testing  of  Mixed  Layer  Models 

Numerous  tests  have  been  completed  or  are  in  the 
planning  stages  for  verification  of  the  forecasting 
abilities  of  the  Garwood,  TOPS,  and  Level-2  Mellor  or  Mellor- 
Durbin  (MD)  models  as  well  as  for  the  other  models  pre¬ 
viously  discussed.  A  review  of  any  mixed  layer  model 
verification  testing  would  not  be  complete  without  at  least 
mention  of  MILE,  the  Mixed  Layer  Experiment,  which  consisted 
of  an  intensive  examination  of  the  upper  ocean  in  the  environs 
of  Ocean  Station  P  during  a  20-day  period  in  autumn  1977. 
Excepting  that  brief  discussion,  the  mainstay  of  the  review 
of  verification  testing  will  encompass  the  three  models 
directly  involved  in  this  study.  (Though  the  MD  model 
is  essentially  the  core  of  the  TOPS  model,  each  will  be 
reviewed  separately,  for  TOPS,  a  first  generation  operational 
model,  has  its  own  inherent  verification  problems.) 

Mixed  Layer  Experiment,  1977:  The  Mixed  Layer  Experiment 
(MILE),  1977,  is  a  noteworthy  example  wherein  data  were 
acquired,  and  compared  initially  against  simulations  produced 
using  the  NI’77  model  and  subsequently  other  model  simulations 
Davis,  et.al.  (1981)  used  velocity  profiles,  measurements  not 
usually  available  in  ML  verification  studies  in  an  attempt 
to  not  only  obtain  a  firm  understanding  of  the  dynamics  and 
thermodynamics  involved,  but  in  so  doing,  to  also  confirm 


the  applicability  of  the  bulk  approach  to  the  study  domain. 
They  failed  to  acquire  such  a  confirmation  in  all  cases 
except  those  involving  strong  storms  wherein  the  upper 
oceanic  shear  became  localized  a  short  distance  below  the 
isothermal  layer.  However,  by  accounting  for  the  vertical 
advection  and  horizontal  pressure  gradient  within  the  domain, 
Davis,  et.al.  (1981)  were  able  to  indicate  a  ’reasonable 
balance'  in  both  the  heat  and  momentum  budgets,  respectively. 
Two  major  difficulties  encountered  were  a)  the  need  to 
acquire  accurate  data  representing  changes  below  the  mixed 
layer  so  as  to  better  understand  the  vertical  distribution 
of  heat  and  momentum  within  the  mixed  layer,  and  b)  the  need 
to  acquire  better  expertise  in  being  able  to  separate  effects 
of  internal  waves,  energetic  inertial  motions  and  internal 
tides  from  responses  of  the  upper  ocean  to  storms.  [Davis, 
et.al.,  1981]  Despite  the  numerous  weaknesses  in  the  study 
described  by  Davis,  et.al.  (1981),  simulations  encouragingly 
seem  to  agree  with  the  observations.  Their  'most  favorable' 
result  (Fig.  4)  was  for  m^  =  0.39  and  m2  =  0.48,  where  m^ 
and  m2  are  tuned  coefficients  determining  'limits  of 
excursion'  and  'rate  of  response',  respectively. 

a.  Review  of  Verification  Testing  of  Three  Models 
The  Level-2  Mellor  model  and  similar  models  have 
been  evaluated  in  numerous  and  varied  simulation  experiments. 
Mellor  and  Durbin  (1975)  obtained  favorable  results  from 
tests  using  five  weeks  of  data  acquired  at  Ocean  Station  Papa 


•  '  t 


48 


MILE  observations  versus  simulations: 

NI  *77,  Dashed  curves:  observed  hour 
average  T(5m)  and  depth  where  T  is  0.1C 
colder  than  T(5m);  solid  curves: 
modeled  values  of  mixed  layer  temperature 
and  depth  [David,  et.al.,  1981] 


(50N,  145W) ,  as  well  as  that  data  acquired  nearby  during 
MILE.  The  model's  ability  to  handle  the  response  of  the 
upper  ocean  to  an  atmospheric  low  pressure  system,  i.e., 
the  observed  versus  the  predicted  SST  during  and  after 
frontal  passage,  was  explored  by  Price  et.al  (1978)  for  the 
Gulf  of  Mexico  region.  The  model's  ability  to  amply  predict 
the  ocean's  response  to  stronger  atmospheric  forcing  was 
verified  in  the  case  of  the  passage  of  hurricane  Eloise  by 
the  fortunate  existence  of  a  NOAA  data  buoy,  EB10,  located 
on  the  track  of  the  hurricane.  Klein  (1980)  used  the  MD 
model  to  simulate  the  variability  of  the  mixed  layer  depths 
in  the  Mediterranean  Sea.  The  model  favorably  simulated 
(Fig.  5)  [Clancy,  et.al.,  1981]  the  diurnal  mixed  layer 
response  observed  during  BOMEX.  (The  Barbados  Oceanographic 
and  Meteorological  Experiment.) 

The  Garwood  model  has  also  undergone,  and  is  still 
undergoing,  verification  testing,  though  by  a  more  limited 
populus .  One  of  the  most  notable  simulations  extended  over 
17  years  of  Ocean  Station  Papa  data.  Not  only  was  this  the 
first  simulation  of  its  kind  in  length,  and  in  the  handling 
of  interannual  variability  in  mixed  layer  response,  but  it 
was  also  the  first  simulation  conducted  for  the  period  of 
the  spring  transition.  (Comparison  of  the  hindcasts  with 
the  observations  lends  encouragement  to  the  ability  of  one¬ 
dimensional  (bulk)  models  to  account  for  a  large  part  of  the 
variance  on  time  scales  from  the  synoptic  to  greater  than  a 
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year.  [Garwood  and  Adamec ,  1982]  Additional  verification 
testing  includes  various  studies  of  model  performance  along 
the  shipping  track  between  San  Francisco  and  Hawaii  [Steiner, 
1981],  and  during  POLEX,  a  component  of  the  North  Pacific 
Experiment.  [Shook,  1980]  Data  assimilation  tests  were 
performed  by  Warrenfeltz  (1980)  and  by  Larsen  (1981). 

Because  the  MD  model  is  at  the  core  of  TOPS,  one 
might  say  that  the  TOPS  model  has  undergone  extensive  testing 
and  evaluation.  Such,  however,  is  not  the  case  in  that 
TOPS  is  a  first  generation  operational  model,  and  as  such  is 
subject  to  inherent  problems  not  experienced  by  other  one¬ 
dimensional  or  quasi-three-dimensional  models.  The  TOPS 
verification  is,  for  the  most  part,  a  test  of  the  quality  of 
the  F NOC  operational  data  base,  the  source  of  not  only  the 
initial  and  upper  boundary  conditions,  but  also  the  verifi¬ 
cation  fields.  These  verification  fields  have  necessarily 
been  extracted  from  the  objective  analysis  or  TOPS  modified 
objective  analysis  scheme,  for  observations  have  rarely  been 
available  near  the  model’s  widely  spaced  gridpoints.  The 
distribution  of  observations  is  neither  uniform  nor 
constant.  The  model  is,  therefore,  subject  to  not  only  a  wide 
range  of  testing  conditions,  but  also  to  a  greatly  fluctua¬ 
ting  source  of  observations.  When  considering  the  advective 
version  of  TOPS,  one  must  also  consider  possible 
three-dimensional  processes. 
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Thus,  meaningful  quantitative  forecast  verification 


seems  reasonable  only  in  regions  containing  an  adequate 
supply  of  data.  The  verification  fields  must  be  obtained  from 
the  same  scheme  used  to  acquire  initial  conditions,  because 
different  schemes  applied  to  the  same  data  set  are  likely  to 
produce  varying  verification  fields.  [Barnett,  et.al.,  1980] 
The  present  FNOC  daily  global  intake  of  data  consists  of 
approximately  200  XBT * s  ,  2,000  'bucket*  temperatures,  and 
approximately  20,000  satellite  SST*s  from  TIROS-N.  This 
apparent  lack  of  sufficient  subsurface  definition  has 
resulted  in  much  of  the  verification  effort  centering  on  the 
SST . 

TOPS  is  an  operational  model,  one  component  of  a 
resource-limited  naval  facility  model  network.  As  such,  it 
neither  has  number  cne  priority  in  the  FNOC  objectives  list, 
nor  can  it  avoid  being  impacted  by  the  ever  evolving  compu¬ 
tational  environment.  The  rate  of  personnel  change,  in 
combination  with  the  shortage  of  oceanographers  on  the  staff, 
has  resulted  in  most  of  the  verification  testing  being  left 
as  a  milestone  to  be  achieved  in  the  future. 

One  of  the  major  changes  which  has  recently 
occurred  at  FNOC,  and  which  will  alter  the  workings  of  TOPS, 
is  the  replacement  of  the  old  atmospheric  primitive  equation 
(PE)  model  with  the  Naval  Operational  Global  Atmospheric 
Prediction  System  (NOGAPS).  This  change  took  place  on 
3  August  1982  midway  during  the  model  simulation  runs  of 
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interest  in  this  thesis.  The  change  was  made  to  improve  the 
higher  priority  atmospheric  products  supplied  to  the  fleet. 

A  study  of  the  effects  of  NOGAPS  on  TOPS  results  is  in  the 
plans  for  the  near  future.  It  is  hoped  that  NOGAPS  will 
provide  more  accurate  heat  fluxes  for  ocean  prediction 
analysis . 

The  analysis  system  has  evolved  from  OTS  (Ocean 
Thermal  Structure)  analysis  system  to  EOTS  (Expanded  Ocean 
Thermal  Structure)  analysis  system  and  finally  to  TEOTS  or 
TOPS-EOTS  (TOPS  modified  EOTS).  The  last  change  in  the 
evolutionary  sequence  was  implemented  in  full  on  22  March  198  3- 
that  is ,  TEOTS  replaced  EOTS  as  a  bottom  boundary  condition 
for  the  atmospheric  models.  Previously,  EOTS  was  the  bottom 
boundary  condition  of  the  atmospheric  models  so  that  there 
was  only  a  one-way  influence  between  the  atmospheric  and 
oceanic  models,  i.e.,  the  atmospheric  model  drove  the 
oceanic  model.  With  the  implementation  of  TEOTS,  the  feed¬ 
back  loop  is  complete.  Thus,  the  bottom  boundary  condition 
of  the  atmospheric  model  is  more  influenced  by  the  physics 
of  the  situation. 

Another  change  to  be  more  fully  discussed  later 
was  a  more  resilient  return  of  the  model  to  climatology.  This 
change  took  place  on  20  September  1982.  [Poliak,  personal 
communication]  Though  not  yet  documented,  the  change 
supposedly  results  in  the  model  forecast  of  SST  remaining 
within  approximately  3C  of  climatology  with  any  deviations 


being  suppressed  within  approximately  a  two-week  period.  (The 
scientific  basis  for  this  criterion  is  unknown.)  (The 
climatology  receives  a  weight  in  the  process  such  that  sharp 
deviations  from  climatology  are  reduced  by  an  increase  of  the 
weighting  applied  to  the  climatological  factor.) 

Most  of  the  published  verification  reports  describe 
TOPS  model  performance  prior  to  these  aforementioned  changes. 
For  example.  Warn  Varnas,  et.al.  (1982)  performed  a  60-day 
simulation  in  the  TRANSPAC  region  of  the  central  North  Pacific 
for  November  and  December  1976.  This  was  done  using  TOPS  in 
conjunction  with  OTS.  Comparison  of  SST  results  from  both 
the  advective  and  nonadvective  versions  of  the  model  revealed 
consistent  skill  over  both  climatology  and  persistence. 

(Fig.  6)  (As  will  be  discussed  later,  results  also  revealed 
a  heat  flux  bias  confirmed  previously  by  comparisons  of  the 
heat  fluxes  with  those  calculated  by  Nate  Clark.)  Clancy, 
et.al.,  (1981)  performed  a  short  time  scale  synoptic  veri¬ 
fication  of  TOPS/EOTS  (nonadvective  version  on  Cyber  175) 
using  both  pattern-of-change  correlation  techniques 
[Dobryshman,  1972],  and  the  apparent  forecast  error 
technique.  The  former  techniques  demonstrated  TOPS’  ability 
to  routinely  produce  real-time  forecasts  of  large  scale 
changes  of  SST.  The  latter  technique  indicated  that,  with 
adequate  data  coverage  and  limited  domains  of  strong 
atmospheric  forcing,  the  model  could  achieve  a  skill  twice 
that  of  persistence. 
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Figure  6.  Skill  of  advective  and  nonadvective  TOPS 
persistence  and  climatology  as  indicated 
TRAMSPAC  data  [Warn-Varnas,  et.al.,  1982] 


The  latest  verifications  of  the  coupled  TOPS 
include  TOPS/TEOTS  vice  TOPS  linked  with  EOTS.  TEOTS  has 
proved  to  be  less  noisy  and  more  consistent  with  the 
predictions  of  TOPS  and  atmospheric  forcing;  this  result 
is  especially  true  in  relation  to  the  spring  transition. 
Composites  of  forecast  verification  statistics  for  May  1981 
indicate  the  skill  of  TOPS  in  forecasting  both  MLD  and  SST 
change  out  to  72  hours.  RMS  errors  indicate  TOPS  betters 
persistence  out  to  (but  not  including)  72  hours  for  MLD, 
but  only  out  to  24  hours  for  SST.  The  latter  deficiency 
is  believed  to  be  the  result  of  the  possible  heat  bias  of 
the  FNOC  primitive  equation  model.  It  is  yet  to  be 
determined  whether  NOGAPS  can  resolve  this  problem. 

C.  STUDY  AREA  AND  ITS  CLIMATOLOGY 

The  domain  investigated  is  (Figs.  1,  2)  located  in  the 
California  Current  System,  just  south  of  the  Mendocino 
Escarpment  (Fig.  7)  in  the  region  from  approximately  37.5  to 
39. 6N  and  125  to  127W.  This  system  is  composed  of  four 
currents:  the  California  Current,  the  Davidson  Current,  the 

California  Undercurrent,  and  the  Southern  California  Current. 
Of  main  influence  in  the  designated  domain  is  the  California 
Current,  the  eastern  boundary  current  of  the  large  anti- 
cyclonic  gyre  centered  near  the  Hawaiian  Islands.  The 
current,  an  extension  of  the  Westwind  Drift  (Sverdrup, 
et.al.,  1942),  occurs  between  the  North  Pacific  atmospheric 
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Study  domain  in  relation  to  bathymetry  off 
Northern  and  Central  California  (Contour 
Interval;  200m) 


high  system  and  the  semipermanent  thermal  low  positioned 
over  central  California.  Geostrophic  and  surface  Ekman 
components  of  the  flow  are  consistent  with  the  temperature 
distribution,  and  the  atmospheric  high  pressure  system  with 
its  associated  wind  stress.  An  anticyclonic  eddy  is  the 
major  feature  depicted  by  the  dynamic  topography  computed 
for  the  CCSI  domain.  The  dynamic  topography  for  CCSII 
clearly  depicts  the  anticyclonic  and  cyclonic  eddies  toward 
the  northeast  and  southwest  regions  of  the  study  domain  which 
dominated  the  region  during  that  time  period  and  into  CCSIII. 
[Reid,  et.al,  1958]  (Figs.  8-10)  Nelson  (1977)  compiled  a 
climatological  description  of  the  surface  wind  stress  and 
wind  stress  curl  over  the  California  Current  System,  by 
computing  the  wind  stress  using  the  bulk  formula  in  one 
degree  squares  from  marine  weather  observations,  primarily 
ship  reports,  archived  by  the  National  Climatic  Center. 

(Figs.  11-13) 

Nelson  encountered  many  of  the  typical  problems  that 
arise  when  dealing  with  multi-source  data  as,  for  example, 
uneven  and  dissimilar  sampling  in  space  and  time.  Only 
approximately  12%  of  the  observations  were  actually  measured. 
The  resulting  values  provide  a  climatological  view  of  the 
area.  The  domain  of  interest  falls  within  a  region  of 
monthly  mean  maximum  wind  stress  values  for  August, 
(approximately  1  dyne /cm* *2)  The  upper  ocean  is  undoubtedly 
influenced  by  such  forcing.  Past  evidence  has  indicated 
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Figure  11.  Distribution  of  observations  per  1-degree 
square;  the  contour  interval  is  2,500 
observations ;  Values  greater  than  5,000 
are  shaded  [Nelson,  1977] 


Figure  12.  Long  term  mean  August  wind  stress 

vectors:  1857-1972.  [Nelson,  1977] 


Figure  13. 


Long  term  me; 
field:  1857- 


that  the  responsiveness  of  the  coastal  regime  to  such 
forcing.  It  is  speculated  that  elements  of  the  coastal 
response  may  be  advected  into  offshore  regions  (and  vice 
versa)  by  the  predominance  of  jets  within  the  region.  (The 
small  negative  value  for  wind  stress  curl  indicates  slight 
open  ocean  downwelling  within  the  study  domain  itself, 
whereas  the  coastal  regions  indicate  strong  upwelling 
tendencies.)  Another  point  of  speculation  is  the  direct 
forcing  of  the  subsurface  eddy  field  by  the  wind  stress  curl 
mechanism  [Mooers,  personal  communication,  1983]. 

Thermodynamical  atmospheric  forcing  fields  are  also 
important  to  mixed  layer  evolution.  Nelson  and  Husby  (1981) 
compiled  monthly  mean  heat  flux  fields  over  the  California 
Current  region.  As  with  the  mean  wind  fields,  the  monthly 
mean  radiative  and  turbulent  heat  fluxes  were  computed  from 
archived  surface  marine  weather  reports  ranging  from 
1921-1972.  The  marine  weather  observations  incorporated  in 
the  heat  flux  calculations  were  of  a  highly  diverse  nature. 
The  distribution  of  the  observations  was  nonrandom  and  highly 
biased  toward  fair  weather.  Also  problems  existed  due  to 
sampling  errors  and  methods  of  computation.  Thus,  large 
errors  were  possible  in  the  heat  flux  estimates.  Errors  of 
10%  in  each  of  the  components  of  heat  exchange  could 
conceivably  have  resulted  in  an  error  of  anywhere  from 
10  to  70%  in  the  net  heat  flux,  Q(N). 

Q(N)=Q(S)  -  Q(B)  -  Q(E)  -  Q(C) 


* 


66 


Despite  the  questionable  magnitude  of  the  long  term  monthly 
means,  however.  Nelson  and  Husby  (1931)  concluded  that  the 
temporal  and  spatial  consistency  of  the  independent  surface 
heat  flux  estimates  ’indicates  that  the  geographic  patterns 
are  realistic  and  significant*.  The  wind  fields  were  of 
similar  or  better  quality. 

The  mean  fields  for  August  of  Q(S),  Q(B),  Q(E),  Q(C), 
and  Q(N)  (Figs.  14-18),  provide  insight  into  heat  gain  and 
loss  within  the  current,  and  more  specifically  within  the 
study  domain.  Values  at  the  domain  center  (38N,  126W),  are: 
Q(S)=216.3  Watts/M**2 
Q(B)=32.77  Watts/M**2 
Q(E)=35.92  Watts /M** 2 
Q(C)=-7.24  Watts/M**2 
Q(N)=154.9  Watts/M** 2 

Standard  errors  of  the  means  evaluated  using 
SE(Q)=SD(Q)/SQRT(N) 

where  SD(Q)  is  the  standard  deviation  and  N  is  the  number  of 
observations  (253  within  the  study  domain  region),  are: 

Q ( S )  : Errors  3 . 80  (Watts /M** 2 ) 

Q(B)  :Error=1.41  (Watts/M**2) 

Q(E)  :Error=3 . 28  (Watts/M**2) 

Q(C)  : Error =1.47  (Watts/M**2) 

Q(N)  : Error  =  5.20  (Watts/M**2) 

As  expected,  the  incoming  solar  incidence  corrected  for 
clouds  and  the  sea  surface  albedo  dominates  the  heat 
budget.  The  effective  back  radiation  and  the  latent  heat 
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Figure  18.  Long  term  mean  August  net  heat 
exchange:  1901-1972  [Nelson 

and  Husby,  1983] 
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term  are  essentially  of  equal  climatological  value  within 
the  region  of  interest.  As  will  be  seen  in  the  following 
discussions,  the  climatological  value  of  net  heat  flux  is 
far  greater  than  that  computed  from  acquired  synoptic 
atmospheric  data  due  to  frequently  overcast  conditions. 

The  climatological  Q(N)  is,  however,  somewhat  less 
(approximately  30  Watts/m**2)  than  the  value  'backed  out* 
from  the  transition  depicted  in  the  acquired  TOPS/TOPS-EOTS 
profiles.  Heat  fluxes  and  wind  stress  values,  the  boundary 
conditions  required  for  the  model  simulations,  are  to  be 
further  discussed  in  Chapter  IV. 

Wind  speed  cubed,  thermocline  strength,  and  mixed  layer 
depth  fields  are  also  included  here  to  demonstrate  the 
spatial  and  seasonal  characteristics  of  wind  generated 
turbulent  kinetic  energy  and  the  associated  upper  ocean  sta¬ 
bility  based  on  Husby  and  Nelson  (1982).  Winds  used  in  the 
computations  of  the  wind  speed  cubed  fields  were  the  same  as 
winds  used  in  the  calculations  of  wind  stress  and  wind 
stress  curl.  The  wind  speed  cubed  fields  were  computed  as  an 
indicator  of  the  turbulent  kinetic  energy  transferred  from 
the  air  to  the  upper  ocean,  a  common  procedure. 

In  producing  the  other  two  fields,  Husby,  et.al.  (1982) 
were  faced  with  some  of  the  typical  problems  one  encounters 
when  working  with  mixed  ''ayers,  more  specifically  with 
interpreting  temperature  profiles  (mainly  XBT).  It  was 
necessary  to  establish  the  criterion  for  the  extraction  of 
the  desired  features.  A  typical  vertical  temperature  profile 
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with  characteristics  of  interest  is  depicted  in  Fig.  3.  (This 
profile  and  some  actual  XBT  profiles  will  be  discussed  in 
more  detail  in  Chapter  III,  Section  A.)  That  BT  profiles 
comprising  their  data  base  were  digitized  at  five  meter 
intervals  plays  an  important  role  in  the  type  of  criterion 
used.  They  chose  to  select  their  mixed  layer  depths  as  the 
upper  limit  of  the  first  five  meter  interval  in  which  the 
temperature  gradient  exceeded  -0.3C/5m  (-0.06C/m).  The 
bottom  of  the  thermocline  was  obtained  by  choosing  the  upper 
depth  of  the  two  successive  five  meter  intervals  within 
which  the  temperature  gradient  did  not  exceed  -0.3 /5m.  The 
strength  of  the  thermocline  was  then  obtained  by  taking  the 
difference  of  the  temperature  at  the  bottom  of  the  thermo¬ 
cline  from  the  mixed  layer  temperature.  These  values  were 
then  summarized  into  one  degree  squares  and  mapped.  (Fig.  19) 
The  practice  of  using  the  wind  speed  cubed  as  an  indica¬ 
tor  as  described  above  does  not  seem  to  apply  to  the  study 
domain  in  August,  for  in  using  an  above  average  estimate  of 
6  m/s  wind,  the  resulting  wind  speed  cubed  is  216  (m/s)**3; 
a  value  quite  a  bit  less  than  that  depicted  in  the  climatology. 
The  mean  observed  mixed  layer  is,  however,  quite  consistent 
with  the  corresponding  climatological  mixed  layer  depicted  in 
Fig.  19.  The  discrepancy  between  turbulence  indicator  (wind) 
and  mixed  layer  depth  leads  one  to  consider  the  probable 
importance  of  convective  activity  and  horizontal  advection 
within  the  region. 
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II.  THE  OBSERVATIONS 


A.  DATA  ACQUISITION 

The  analyzed  data  were  acquired  aboard  the  R/V  ACANIA 
in  support  of  the  Pilot  Ocean  Prediction  Study  of 
California  Current  Eddies,  which  has  since  become  known  as 
the  OPTOMA  program  sponsored  by  ONR.  The  cruises  designated 
CCS1  and  CCS2  (Legs  I  and  II)  covered  time  intervals  8  to 
12  March,  30  July  to  5  August,  and  9  to  14  August  1982, 
respectively.  (For  the  convenience  of  this  thesis,  CCS1 
will  be  called  CCSI,  CCS2  legs  I  and  II  will  be  called  CCSII 
and  CCSIII,  respectively.)  Each  had  as  its  domain  of 
interest,  a  region  off  Point  Arena  centered  at  approximately 
38N ,  126W.  CCSII  was  devoted  to  sampling  the  largest  domain, 
approximately  180  km  square  with  line  spacing  of  ca.  44km. 
and  station  spacing  of  approximately  a  quarter  of  the 
baroclinic  Rossby  radius  of  deformation,  i.e. ,9  km.  This 
allowed  for  better  resolution  over  a  scale  neither  so  large 
as  that  of  the  CalCOFI  data  nor  as  small  as  that  used  during 
CCSI.  (120  km  with  line  spacing  of  ca.  18  km  and  station 
spacing  of  9  km.)  CCSIII  was  for  the  most  part  a  resampling 
of  the  CCSII  grid.  Station  spacing  on  a  whole  was  maintained 
at  small  enough  intervals  (ca.  9  km  for  August  cruises)  to 
minimize  aliasing  by  the  smaller  scale  activity. 


On  all  three  cruises,  continuous  sampling  was  achieved 
using  the  shipboard  DAS  (Digital  Data  Acquisition  System) . 

T-4  XBT's  (expendable  bathythermographs)  maximum  depth  of 
450  m,  were  launched  at  each  station.  (Figs.  1,  2)  189  XBT's 

were  deployed  on  CCSI,  while  approximately  twice  that  number, 
353,  were  deployed  in  the  combined  latter  two  cruises.  No 
CTD  casts  were  made  during  CCSI.  Interspersed  along  the 
track  on  CCSII  and  CCSIII  were  500  and  1500m  CTD  casts. 

During  CCSII  and  CCSIII,  100  CTD  casts  were  made  prior  to 
equipment  failure  due  to  rough  weather.  In  addition  to 
logging  these  vertical  profiles,  DAS  stored  station  and 
interstation  "underway"  values  of  position,  time,  ship's 
speed,  ship's  roll,  and  surface  atmospheric  data;  such  as, 
dewpoint  temperature,  air  temperature,  barometric  pressure, 
wind  direction,  and  wind  speed.  Fractional  cloud  cover  was 
recorded  in  the  ship's  log  on  an  hourly  basis  by  the  deck 
watchstander . 

Winds  for  each  day  of  the  August  cruises  in  the  vicinity 
of  the  cruise  tracks  (nine  gridpoints)  were  supplied  by  FNOC. 
(Figs.  20,  21)  In  addition,  TOPS  and  TOPS-EOTS  produced 
ocean  temperature  profiles  associated  with  eight  of  the  nine 
points  of  interest  which  were  archived  by  FNOC  for  later 
evaluation  with  in  situ  data.  (The  ninth  point  was  located 
on  land.) 

Satellite  IR  imagery  was  available  from  NOAA/NESS 
(Redwood  City)  for  cruise  planning  and  preliminary  analysis. 
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Higher  quality  imagery  was  supplied  by  SOSF  (Scripps  Ocean 
Satellite  Facility)  for  post  analysis. 


B.  DESCRIPTION  OF  DATA  AVAILABLE 

1.  CCS I 

The  XBT  data  are  of  primary  interest  from  the  March 
cruise.  189  XBT's  were  deployed,  of  which  180  were 
'successful’.  Each  trace  was  digitized  by  FNOC  because 
DAS  had  malfunctioned.  Final  selection  left  105  XBT's 
within  the  grid  region.  Data  files  consisted  of  date,  time, 
latitude,  longitude,  and  vertical  profiles  of  depth  versus 
the  corresponding  temperature  and  climatologically-based 
salinity. 

2.  CCSII  and  CCSIII 

These  cruises  contained  more  data  of  present  interest. 
The  XBT  count  was  257  after  editing.  Digitization  by  FNOC 
in  this  case  was  not  necessary.  Of  these  XBT's,  198  were 
located  within  the  grid. 

The  data,  similar  to  that  of  CCSI  with  the  exception 
that  the  salinity  was  not  climatogically-based ,  were  recorded 
on  DAS  cartridges.  Data  from  these  tapes  were  subsequently 
made  accessible  from  CMS  (Coversational  Monitor  System)  on 
the  NPS  IBM  3033.  The  original  traces  and  XBT  logs  were  used 
for  comparison  and  error  checks. 

DAS  log  data,  stored  and  accessed  in  the  same  way  as 
the  XBT  data,  though  in  the  standard  DAS  header  file  format. 
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were  the  source  of  the  majority  of  the  atmospheric  data 
acquired  during  CCSII  and  CCSIII.  Numerous  runs  in  various 
editing  modes  were  made  to  delete  blocks  of  zeros,  erroneous 
values,  and  data  not  applicable  to  the  grid  region,  (see 
Appendix  A) 

Additional  atmospheric  data,  such  as  cloud  cover, 
were  recorded  on  an  hourly  basis  in  the  ship's  log.  Any 
additional  data  pertaining  to  the  ocean,  such  as  bucket 
temperatures,  were  logged  in  the  CTD  log  or  on  the  XBT 
traces . 

Finally,  atmospheric  as  well  as  oceanic  analyses  in 
the  vicinity  of  the  experimental  domain  were  archived  by  FNOC. 
They  consisted  of  daily  wind  values,  one  for  each  of  the 
nine  selected  FNOC  gridpoints,  as  well  as  daily  analyses  and 
24-hour  forecasts  of  the  ocean  temperature  profiles  for  the 
eight  gridpoints  located  over  the  ocean.  Since  TOPS/TOPS-EOTS 
runs  on  a  Northern  Hemispheric  polar  stereograDhic  grid  with 
a  grid  spacing  of  approximately  300  km  in  the  midlatitudes, 
it  was  not  possible  to  obtain  analyses  from  numerous  points 
in  the  cruise  region.  In  fact,  of  the  eight  points  located 
over  the  ocean,  none  of  them  fell  directly  within  the  cruise 
domain.  (Figs.  20,  21)  The  use  of  the  standard  FNOC  gridp~int 
was  an  unnecessary  limitation  that  was  imposed  for  convenience. 
(Subsequently,  it  was  learned  that  FNOC  has  the  capacity  to 
interpolate  their  analyses  to  produce  a  vertical  temperature 
profile  at  a  specified  point.)  The  analyses  and  TOPS 
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predictions  used  in  this  study  were  from  the  point  closest 
to  the  study  domain.  Based  on  the  SST  climatology  (Fig.  22), 
the  chosen  FNOC  gridpoint  at  39.Q5N,  127. 75W  was  not 
expected  to  have  a  SST  much  greater  (at  most  1C)  than  that 
of  the  XBT’s  within  the  cruise  area. 

C.  FORCING  FUNCTIONS  FOR  THE  GARWOOD  AND  MELLOR  LEVEL-2.5 

MODELS 

The  version  of  the  Garwood  model  utilized  required  an 
initial  vertical  temperature  profile  as  well  as  three-hourly 
values  of  air  temperature,  dewpoint  temperature,  wind  speed, 
and  cloud  cover  for  the  computation  of  the  boundary  conditions 
to  be  used  in  both  the  Garwood  and  the  Mellor  2.5  models: 
the  momentum  and  heat  fluxes. 

For  the  first  run,  initialization  was  done  with  the  FNOC 
analyzed  ocean  temperature  profile  at  39.  ON,  127. 75W.  This 
initialization  permitted  comparison  of  the  forecasted  TOPS/ 
TOPS-EOTS  ocean  temperature  profile  with  that  obtained  from 
the  Garwood  model  over  the  same  time  span.  All  other  runs 
were  made  using  XBT  96  as  the  initial  temperature  profile. 

The  three-hourly  air  temperature  was  obtained  by 
averaging  the  DAS  underway  data.  Averages  were  made  over  6 
to  299  values.  The  resulting  standard  deviations  range  from 
a  high  of  1.67C  to  a  low  of  0.0005C.  Only  three  values  had  a 
standard  deviation  over  1.0C.  Averaging  was  done  under  the 
assumption  that  the  air  temperature  varied  little  over  the 
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Figure  22.  Average  SST  for  last  half  of  July  1982 


relatively  small  cruise  grid  at  any  given  time.  In  addition, 
it  was  assumed  in  the  case  of  missing  or  erroneous  data 
that,  at  a  given  time  of  day,  the  values  did  not  vary  much 
over  the  course  of  a  couple  of  days.  Thus,  for  missing 
0001-0300  readings,  the  0300  input  to  the  model  was 
assumed  equal  to  that  of  the  previous  day.  (Fig.  23) 

The  dewpoint  temperatures  were  to  be  acquired  in  a 
similar  manner,  however,  faulty  instrumentation  precluded 
this.  Dewpoint  temperatures  were  instead  computed  using  a 
technique  described  by  Bolton  (1980).  (Fig.  24)  To  use  the 
method,  relative  humidities  (also  unavailable)  were  required. 

As  recommended  by  Davidson  (personal  communication),  a 
relative  humidity  of  90%  was  assumed  throughout  to  correspond 
with  the  foggy/hazy  conditions  which  existed  intermittently 
through  the  cruises. 

Wind  speeds  (Fig.  25)  from  the  nearest  FNOC  point  were 
used  consistent  with  the  location  of  initial  values.  The 
region  of  interest  was  considered  to  be  too  small  to  require 
interpolation  of  winds.  (Fig.  26)  Ship  winds  available  from 
DAS  are  relative  winds,  but  true  winds  could  be  calculated 
by  correcting  for  ship's  velocity.  This  was  not  necessary 
for  the  runs  described.  Instead,  the  applicability  of  the 
FNOC  winds  was  checked  by  comparing  them  with  the  DAS  winds 
recorded  at  CTD  stations  because  the  ship  was  then  'stationary', 
(s'ee  Appendix  A) 
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Time  series  of  three  hourly  averaged  DAS  recorded 
temperatures  used  in  computing  forcing  functions  f 
the  Garwood  and  Mellor  Level-2.5  model  simulations 


Time  series  of  calculated  dewpoint  temperatures 
used  to  compute  forcing  functions  for  the  Garwood 
and  Mellor  Level-2.5  model  simulations 


eries  of  FNOC  analyzed  wind  components  and 
for  39.05N,  127. 72W 


Sea  surface  temperatures  were  obtained  from  individual 
XBT's  along  the  track.  (Fig  27)  The  XBT's  used  were  those 
launched  within  half  an  hour  of  the  necessary  input  time. 

The  final  input,  cloud  cover  in  octals,  was  acquired 
from  hourly  observations,  and  averaged  over  three-hour 
intervals.  (Fig.  28) 

The  boundary  conditions,  i.e.,  heat  and  momentum  fluxes, 
applied  were  computed  within  the  program.  (Figs.  29,  30)  As 
will  be  discussed  later,  the  average  daily  net  heat  flux 
forcing  the  Garwood  and  Mellor  models  over  the  thirteen  days 
was  +4.25  W/m**2.  (where  positive  is  upward  flux,  i.e.,  the 
ocean  was  losing  heat  to  the  atmosphere) 

Additional  information  required  by  the  Garwood  model 
includes  various  coefficients;  such  as,  the  extinction 
coefficient  and  the  drag  coefficient.  For  the  most  part,  they 
were  set  from  previous  model  usage.  The  Garwood  model  con¬ 
tains  two  'tuning  parameters'  (resulting  in  essentially  two 
degrees  of  freedom)  which  had  to  be  addressed  prior  to  model 
simulations  with  each  'new  host  point'  [Garwood,  et.al., 

1982]  The  parameters,  AM  and  AZR  had  already  been  set  upon 
receipt  of  the  model.  (A  tuning  analysis  is  presently 
underway  within  the  Ocean  Prediction  Laboratory.  The  procedure 
and  results  will  probably  be  documented  at  some  future  time.) 

Coefficients,  such  as  the  diffusion  coefficients,  used 
in  the  Mellor  level- 2. 5  model  are  'universal'  constants,  and 
so  required  no  alterations. 


Figure  27.  Time  series  of  XBT  S*3T  used  to  compute  the  forcing 
functions  for  the  Garwood  and  Mellor  Level-2.5  model 


III.  STATISTICS 


A.  DETERMINATION  OF  THE  OBSERVED  VARIABLE  MIXED  LAYER  DEPTH 

Obtaining  mixed  layer  depths  from  DAS  or  FNOC  digitized 
profiles  or  even  from  the  original  XBT  traces  is  no  easy 
task.  Criteria  employed  vary  greatly,  and  there  is  a  ten¬ 
dency  to  change  one’s  criterion  with  each  new  exception  that 
arises.  Further  complexity  occurs  with  the  desire  to  isolate 
new  (shallow)  mixed  layer  depths  in  addition  to  the  numerous 
old  (deep)  ones. 

The  new  mixed  layer  is  a  layer  formed  over  the  old  one 
due  to  the  mixing  being  insufficient  to  reach  the  level  of 
the  old  (previously  formed)  mixed  layer.  This  may  be  the 
result  of  weak  winds  and  surface  heating.  The  old  mixed 
layer  will  remain  in  place  until  turbulent  diffusion  of  the 
interface  or  some  other  process  dissipates  it  or  strong  mixing 
deepens  the  layer.  The  main  reason  for  this  is  that  the 
Second  Law  of  Thermodynamics  does  not  permit  the  unmixing  of 
the  fluid.  An  ’excellent*  example  of  an  old  mixed  layer 
recorded  during  CCSI  is  depicted  in  Fig.  31.  (This  layer  is 
also  the  new  mixed  layer  for  it  has  replaced  the  old  one 
which  would  have  existed  at  an  earlier  time.)  An  old  mixed 
layer  extending  to  approximately  63m  with  an  overlying  new 
mixed  layer  extending  to  approximately  41m  is  depicted  in 
Fig.  32.  In  Fig.  3,  h^  would  be  considered  a  new  mixed  layer 
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overlying  the  old  winter  mixed  layer,  hj.  New  mixed  layers 
do  not  depend  solely  on  seasonal  or  synoptic  changes. 

Numerous  cases  of  diurnal  mixed  layers  were  noted  wherein 
the  cooling  of  the  surface  at  night  caused  the  upper  water 
column  to  become  unstable  and  consequently,  mixed,  forming  a 
new  mixed  layer.  Meanwhile,  during  the  day,  the  surface  is 
heated,  possibly  producing  a  warm  shallow  layer  over  the  old 
mixed  layer,  i.e.,  a  new  thermocline  is  formed. 

New  mixed  layers  may  also  be  formed  under  a  large  number 
of  other  conditions.  A  storm,  for  example,  will  most  likely 
cause  a  new  mixed  layer  to  form  which  will  quite  possibly 
be  deeper  than  the  old  one.  In  this  case,  the  new  mixed 
layer  supplants  the  old  one,  and  becomes  the  old  mixed  layer, 
too . 

Eventhough  the  old  mixed  layers  are  of  importance  to  Navy 
applications,  for  their  role  in  surface  acoustic  ducting,  the 
new  mixed  layers  are  the  layers  of  importance  in  mixed  layer 
modeling,  i.e.,  that  is  the  layer  predicted  and  modified  by 
some  bulk  (layered)  models  in  use  today.  Thus,  one  must 
find  a  suitable  method  for  defining  both  of  these  significant 
features,  as  well  as  sea  surface  temperature,  depicted  in  an 
XBT  profile. 

The  first  attempt  to  obtain  mixed  layer  depths  from  the 
temperature  depth  profiles  was  made  by  computer  comparison 
of  three  sequential  points  on  the  trace  for  a  0.2C  change. 
(The  value  0.2C  was  chosen  for  that  is  the  accuracy  of  the 


new  XBTs ;  more  precisely,  0.12C-0.31C  [Heinmiller,  1982], 

It  is  the  value  presently  in  use  in  producing  the  TOPS/ 
TOPS-EOTS  contoured  mixed  layer  depth  product.)  This 
technique  proved  to  be  satisfactory  except  in  cases  of  a 
surface  thermocline,  a  very  weak  new  mixed  layer,  no  mixed 
layer,  and  an  old  mixed  layer  with  a  noticeable  gradient 
(though  still  possessing  the  characteristic  'knee').  A 
similar  technique  was  used  by  Thompson  (1975)  with  the 
additional  technique  of  starting  his  sequential  comparison  at 
three  meters  depth  to  avoid  the  surface  thermocline  problem. 
With  either  procedure  one  would  expect  to  obtain  a  depth,  no 
matter  whether  a  mixed  layer  existed  or  not.  (The  depth  for 
the  0.2C  change  obtained  may  be  in  the  middle  of  a  gradient.) 
Such  was  not  the  case,  however,  due  to  the  digitization 
scheme  used  by  FNOC  wherein  only  depths  of  inflection  points 
were  listed.  Comparison  of  the  surface  temperature  with  the 
temperature  of  the  first  inflection  point,  the  point  at  the 
base  of  the  shallow  surface  thermocline,  may  yield  a  greater 
than  0.2 C  change.  This  would  in  turn  lead  to  0.0m  being 
denoted  the  mixed  layer  depth  when  a  very  noticeable  mixed 
layer  might  be  apparent  below  the  shallow  thermocline  or 
possibly  due  to  the  response  time  of  the  XBT  recorder. 

These  possible  deficiencies  made  it  necessary  to  return 
to  the  original  traces  and  logs.  Each  trace  was  examined  and 
redefined  in  terms  of  temperatures  and  depths  of  the  features 
of  interest,  for  example,  surface  temperature,  the  temperature 
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and  depth  of  the  base  of  the  surface  thermocline,  the 
temperature  and  depth  of  each  sequential  possible  mixed 
layer,  and  the  temperature  and  depth  of  a  deeper  point  chosen 
to  define  the  seasonal  thermocline.  Final  mixed  layer  depth 
values,  as  well  as  other  values,  were  chosen  by  comparison 
of  the  computerized  analysis,  the  traces,  and  the  XBT  logs. 
The  computerized  analysis  was  given  the  most  weight,  followed 
by  results  obtained  from  inspection  of  the  traces.  The  XBT 
logs  were  given  the  least  emphasis  due  to  their  variation 
with  the  experience,  knowledge,  and  judgment  of  the  observer 
and  the  not  always  optimal  conditions  under  which  they  were 
recorded. 

B.  STATISTICAL  CHARACTERIZATION 

The  several  variables  derived  from  the  XBT  profiles  by 
the  techniques  previously  mentioned  were  used  as  a  basis  for 
a  statistical  characterization  of  each  cruise.  Basic  pro¬ 
perties  such  as  the  mean,  minimum,  maximum,  median,  and 
standard  deviation  were  examined.  They  are  summarized  in 
Table  III,  Two  types  of  correlations  coefficients  were 
calculated,  as  well  as  Student-t  confidence  intervals, 
regression  equations,  and  histograms  and  maps  of  SST,  bucket 
temperature,  old  mixed  layer  depth  (OMLD),  old  mixed  layer 
depth  temperature  (OMLDT),  new  mixed  layer  depth  (NMLD)  and 
new  mixed  layer  depth  temperature  (NMLDT).  Only  a  limited 
selection  of  statistical  quantities  will  be  presented  to 
convey  the  character  of  the  upper  ocean  during  the  various 
cruises. 
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The  first  variables  to  be  discussed  are  the  SST's  for 
the  three  cruises  and  the  bucket  temperatures  from  the  first 
cruise.  A  total  of  105  XBT's  were  available  from  the  CCSI 
grid  for  the  statistical  characterization.  Both  XBT  SST  and 
bucket  temperature  were  recorded  with  a  comparison  revealing 
the  XBT  SST  field  was  overall  slightly  warmer  than  the  bucket 
temperature  field,  i.e.,  0.4C.  (Fig.  33,  Table  III)  A 
regression  equation  or  least  squares  linear  equation  of  the 
bucket  temperature  with  respect  to  the  XBT  SST  was  calculated 
to  be:  Y= . 678  +  .926(X1)  where  X1=XBT  SST.  The  R  squared 
value,  a  measure  of  how  well  the  regression  equation  fits 
the  data,  was  79.6%.  (100%  is  a  perfect  fit.  The  R  squared 
value  is  defined  by  (100) (sum  of  squares  due  to  regression/ 
total  sum  of  squares.)) 

The  two  correlation  coefficients  computed  were  the  Pearson 
product  moment  coefficient  and  the  autocorrelation  function 
with  respect  to  distance.  The  former  calculated  for  XBT  SST's 
versus  bucket  temperatures,  resulted  in  a  large  correlation 
coefficient,  .892.,  though  not  'large  enough'  considering 
the  fact  that  the  values  should  be  exactly  the  same  (except 
for  measurement  noise).  The  spatial  autocorrelation  function 
was  computed  using  binned  XBT  pairs.  Out  of  210  km,  the 
radial  resolution  for  XBT  pairs  within  a  bin  was  15  km. 

Beyond  210  km,  it  was  increased  to  25  km.  XBT  SST  for  CCSI 
(Fig.  34)  had  a  correlation  distance  of  approximately  30  km, 
using  a  0.2  correlation  as  a  reference.  The  correlation 
distance  for  the  bucket  temperature  was  slightly  less,  i.e., 

27  km.  (Fig.  35) 
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Figure  33.  Histograms  of  CCSI  bucket  and  SST, 
CCSII  SST,  and  CCSIII  SST 
(Temperatures  in  degrees  Celsius) 
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CCS1  XBT  SST  RUTOCORRELRTION  FUNCTION 


SHIYd  dO  'ON 


35.  Autocorrelation  function  for  CCSI  XBT  SST  with  respect  to 
distance:  results  depicted  are  for  detrended  data  and 

’raw’  data  (with  trend);  also  indicated  are  the  number  of 
data  pairs 


TABLE  III 


BASIC  STATISTICS 


CCSI 

CCSII 

CCSIII 

SST 

Bkt 

SST 

SST 

Mean 

11.5 

11.3 

15.9 

16.3 

Stddev 

.302 

.314 

.91 

.703 

Min 

10.9 

10.4 

12.9 

13.9 

Max 

12.4 

12.0 

18.0 

17.1 

Median 

11.5 

11.3 

15.8 

16.5 

T  90% 

(11.4,11.5)  (11.3,11. 

4)  (15.7,16.1) 

(16.1,16.4) 

C.I.* 

NMLDT 

NMLDT 

NMLDT 

Mean 

11.3 

15.4 

15.9 

Stddev 

.324 

.822 

.7  44 

Min 

10.8 

12.3 

11.9 

Max 

12.2 

16.7 

16.9 

Median 

11.3 

15.6 

16.2 

T  90% 

(11.3,11 

.4) 

(15.3,15.6) 

(15.8,16.1) 

C.I.* 

NMLD 

NMLD 

NMLD 

Mean 

39.1 

15.1 

16.3 

Stddev 

25.5 

8.3 

8.1 

Min 

1.0 

4.0 

5.0 

Max 

82.0 

55.0 

50.0 

Median 

33.0 

13.0 

15.0 

T  90% 

(34.9,43 

.4) 

(13.6,16.6) 
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OMLDT 
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OMLDT 

Mean 

11.2 

14.7 

15.5 

Stddev 

.287 
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.801 
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10.5 
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11.2 

14.8 

15.7 

T  90% 

(11.1,11 

.2) 

(14.6,14.9) 

(15.4,15.6) 

C.I.* 

OMLD 

OMLD 

OMLD 

Mean 

64.8 

37.8 

29.5 

Stddev 

12.0 

13.4 

12.2 

Min 

27.0 

6.0 

8.0 

Max 

95.0 

69.0 

65.0 

Median 

63.5 

38.0 

28.0 

T  90% 
C.I.* 

(62.8,66 

.7) 

(35.5,40.1) 

(27.5,31.6) 

Unites 

in  table:  all 

temperatures:  degrees  Celsius 

all 

depths ; 

meters 

T  Interval  With  90%  Confidence 


*  T  90%  C.I. : 


During  the  latter  two  cruises,  bucket  temperatures  were 
not  recorded  routinely.  The  CCSII  data  set  consisted  of  97 
XBT’s.  As  was  to  be  expected,  the  SST  values  were  much 
higher  than  those  recorded  for  the  March  cruise.  In  fact, 
the  minimum  CCSII  temperature,  12. 9C,  exceeded  the  maximum 
CCSI,  12. 4C.  The  range  of  temperatures  for  CCSII,  and  thus 
the  standard  deviation  in  this  case  was  quite  large.  A 
reason  for  the  significant  range,  12.9  -  18. OC,  was  the 
presence  of  cold  eddies  and  cool  anomalies  near  the  peak 
of  summertime  warming. 

The  autocorrelation  functions  for  SST  for  CCSII  were 
plotted  with  and  without  a  linear  trend  removed.  Fig.  36. 

This  was  done  to  get  some  idea  of  the  existence  and  effects 
of  a  trend  in  the  data.  Also  shown  is  a  95%  confidence 
region  computed  by  +/-  sqrt(2)/N  where  N  is  the  number  of 
paris  in  the  bin.  (Using  this  N  instead  of  the  total  N 
results  in  a  more  conservative  confidence  region.)  Results 
indicate  a  slight  increase  of  the  correlation  distance 
compared  to  the  March  cruise,  i.e.,  32  km  for  the  detrended 
data  and  38  km  without  the  trend  removed. 

The  mean  SST  for  CSSIII  was  still  higher,  16. 2C,  as  might 
be  expected  for  later  in  August.  The  range,  however,  fit 
within  that  of  the  second  cruise.  This  may  be  related  to 
the  fact  that  CCSIII  did  not  cover  as  large  a  region  as  CCSII 
The  correlation  distance  for  this  cruise  was  40  km  for  the 
detrended  data  and  43  km  without  the  trend  removed.  (Fig.  37) 
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re  36.  Autocorrelation  function  For  CCSII  XBT  SST 

with  respect  to  distance;  Results  depicted  are 
for  detrended  data  and  'raw'  data;  Also 
indicated  are  the  number  of  data  pairs  and  the 
95%  confidence  limits 
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Autocorrelation  function  for  CCSIII  XBT  SST . 
with  respect  to  distance;  Results  depicted  are 
for  detrended  data  and  ’raw*  data;  Also 
indicated  are  the  number  of  data  pairs  and  the 
95%  confidence  limits 


This  might  indicate  a  more  homogeneous  area  as  compared  to 
the  larger  grid  transited  during  CCSII. 

The  new  mixed  layers  for  the  March  cruise  ranged  from 
0.0  m  (when  no  new  mixed  layer  was  observed)  to  82  m.  (Table 
III)  The  largest  Pearson  products  moment  correlation  obtained 
was  for  the  NMLD's  and  OMLD's.  This  resulted  due  to  a  number 
of  cases  when  the  NMLDwas  identical  to  the  OMLD,  i.e.,  in 
cases  of  intense  turbulence. 

For  CCSII,  84  XBT's  possessed  a  NMLD.  The  NMLD's  were, 
however,  quite  small,  approximately  half  of  those  noted  in 
CCSI,  (Fig.  38).  CCSIII  had  94  instances  of  NMLD's.  The 
distribution  in  the  histogram  appears  to  be  Raleigh  in 
nature . 

Histograms  of  the  OMLD's  for  the  cruises  are  shown  in 
Fig.  39.  Most  noticeable  is  the  shift  toward  shallower 
mixed  layers  with  time.  In  fact,  the  mean  depth  for  the  March 
cruise  is  a  factor  of  two  greater  than  the  CCSIII  mean  OMLD. 
The  Pearson  product  moment  correlation  coefficients  for  the 
OMLD  in  association  with  the  majority  of  the  other  variables 
showed  no  dominant  results.  The  largest  coefficient  obtained 
was  with  respect  to  the  NMLD;  the  probable  reason  for  this 
was  already  discussed.  For  CCSI,  the  correlation  distance 
was  merely  15  km.  (Fig.  40)  For  CCSII,  both  the  data  with 
and  without  the  trend  removed,  have  a  correlation  distance 
of  only  approximately  14  km.  (Fig.  41)  CCSIII  possessed  old 
mixed  layers  that  were  more  highly  correlated  than  those 
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Figure  38. 


Histograms  of  NMLD  (meters)  for  CCSI , 
CCSII,  and  CCSIII. 
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Figure  40.  Autocorrelation  function  for  CCSI  old  mixed 
layer  depth 


Figure  41.  Autocorrelation  function  for  CCSII  old  mixed 
layer  depth 


w 


previously  examined.  (Fig.  42)  The  correlation  distance  for 
the  detrended  data  was  approximately  35  km;  without  the 
trend  removed  the  correlation  distance  was  approximately  50 
km.  The  increase  in  the  correlation  distance  from  CCSII  to 
CCSIII  may  be  due  to  the  fact  that  the  former  cruise  covered 
a  wider  and  more  diversified  domain  than  did  the  latter. 

Contours  of  OMLD  are  included  for  CCSII  and  CCSIII  to 
give  an  idea  of  the  variability  of  mixed  layers  (a  quasi¬ 
synoptic  depiction),  and  to  bring  to  light  a  common  problem 
encountered  in  trying  to  contour  'real'  data.  The  data  must 
be  placed  on  a  regularly  spaced  grid  prior  to  using  contouring 
graphics  packages.  One  would  think  the  answer  to  the  problem 
would  be  to  merely  interpolate  the  data  to  the  new  grid. 
Problems  arise,  however,  from  data  set  to  data  set  depending 
on  the  noise  or  gradients  in  the  data.  Certain  schemes  seem 
to  work  well  for  some  sets  and  not  to  work  for  others.  The 
enclosed  computer  plotted  contoured  fields  are  the  result  of 
interpolation  by  trying  to  fit  a  five  degree  polynomial  to  the 
data.  (Fig.  43)  The  results  do  not  have  much  worth  as  seen 
when  compared  with  the  hand  contours.  (Fig.  44)  Only  the 
locations  of  some  of  the  minimum/maximum  centers  are  retained. 
In  the  future,  the  technique  of  optimal  interpolation  should 
be  applied  to  the  data  for  improved  interpolated  values  and 
consequently,  more  appropriately  smoothed  contoured  fields. 
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Figure  42.  Autocorrelation  function  for  CCSIII  old  mixed 
layer  depth 
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Figure  43 


Computer  Contoured  old  mixed  layer 
depths  from  CCSIII  data  were 
interpolated  to  a  regular  rotated 
grid  and  contoured;  rotation  origin: 
37.33N,  126. 25W;  subplot  origin: 

38. ION,  126. 38W 
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IV.  ANALYSIS  OF  OBSERVED  AND  SIMULATED 

- TTHmiTingTSTRUCTTOE - 

The  first  run  was  of  the  Garwood  model  initialized  with 
the  TOPS  predicted  profile  for  1  August.  The  initial  mixed 
layer  depth  was  set  at  20m,  and  the  initial  SST  chosen  was 
17. 7C  (the  SST  of  the  TOPS  analysis  profile).  The  24  daily 
profiles  were  averaged  so  as  to  produce  one  'representative* 
profile  per  day.  The  average  profiles  along  with  their 
standard  deviations  compose  Figs.  45-58.  The  standard 
deviation  was  nonzero  only  in  the  region  at  and/or  above  the 
base  of  the  mixed  layer,  because  the  Garwood  model  does  not 
alter  regions  below  the  furthest  extent  of  the  mixed  layer 
base.  For  the  13-day  cruise  period,  the  Garwood  model 
cooled  the  mixed  layer  by  approximately  two  thirds  of  a 
degree,  and  deepened  the  mixed  layer  to  approximately  37m. 

A  substantial  amount  of  entrainment  occurred  from  2  to  3 
August,  with  the  simultaneous  increase  in  wind  speed  (7.8  m/s), 
and  decrease  in  the  SST  at  00Z.  Thus,  there  was  increased 
forced  convective  activity  due  to  increased  wind  stirring  and 
increased  free  convection  due  to  surface  ccoling.  Significant 
mixing  also  occurred  over  the  time  spans:  7  to  8,  9  to  10, 

12  to  13  August  with  some  entrainment  also  occurring  from 
11  to  12  and  10  to  11  August,  but  not  as  pronounced  as  that 
during  the  previously  stated  times.  Not  much  occurred  during 
the  other  time  spans. 
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AYE  MODEL  PROFILE,  AUG  1 


TEttPERRTURE 

0.0  5.0  10.0  15.0  20.0  2S.0 


Average  1  August  profile  from  Garwood 
fnodel  .  (Temperature  in  degrees  Celsius 
Initialized  with  FNOC  analyzed  profile 


Figure  47.  Average  2  August  profile  from  Garwood 
model,  Initialized  with  FNOC  analyzed 
profile  (Temperature  in  degrees  Celsius) 
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AVE  MODEL  PROFILE,  AUG  3 
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8.  Average  3  August  profile  from  Garwood 
model  Initialized  with  FN0C  analyzed 
profile  (Temperature  in  degrees  Celsius) 
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'e  49.  Average  4  August  profile  from  Garwood 
model.  Initialized  with  FNOC  analyzed 
profile  (Temperature  in  degrees  Celsius) 
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AVE  MODEL  PROFILE,  AUG  5 


TEMPERATURE 
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Figure  50.  Average  5  August  profile  from  Garwood 
model.  Initialized  with  FNOC  analyzed 
profile  (Temperature  in  degrees  Celsius) 
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AYE  MODEL  PROFILE,  AUG  6 


TEMPERATURE 


Figure  51.  Average  6  August  profile  from  Garwood 
model.  Initialized  with  FNOC  analyzed 
profile  (Temperature  in  degrees  Celsius) 


124 


AYE  MODEL  PROFILE,  AUG  7 


TEMPERATURE 


Average  7  August  profile  from  Garwood 
model.  Initialized  with  FNOC  analyzed 
profile  (Temperature  in  degrees  Celsius) 


AYE  MODEL  PROFILE,  AUG  8 


temperature: 


Figure  53.  Average  8  August  profile  from  Garwood 
model.  Initialized  with  FNOC  analyzed 
profile  (Temperature  in  degrees  Celsius) 
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AVE  MODEL  PROFILE,  AUG  9 


TEMPERATURE 
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Figure  54.  Average  9  August  profile  from  Garwood 
model  Initialized  with  FN0C  analyzed 
profile  (Temperature  in  degrees  Celsius) 
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AYE  MODEL  PROFILE,  AUG  10 


TEMPERATURE 

0.0  5.0  10.0  1S.0  20.0  2S.0 


Figure  55.  Average  10  August  profile  from  Garwood 
model  Initialized  with  FN0C  analyzed 
profile  (Temperature  in  degrees  Celsius) 
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Figure  56.  Average  11  August  profile  from  Garwood 
model  Initialized  with  FNOC  analyzed 
profile  (Temperature  in  degrees  Celsius) 
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AVE  MODEL  PROFILE,  AUG  12 


TEMPERATURE 


Figure  57.  Average  12  August  profile  from  Garwood 
model  Initialized  with  FNOC  analyzed 
profile  (Temperature  in  degrees  Celsius) 
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AVE  MODEL  PROFILE,  AUG  13 


.gure  58.  Average  13  August  profile  from  Garwood 
model.  Initialized  with  FNOC  analyzed 
profile  (Temperature  in  degrees  Celsius) 
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The  time-depth  contours  of  temperature  in  Fig.  59  are 
another  way  of  indicating  the  deepening  and  shoaling  of  the 
mixed  layer  depth.  It  is  quite  obvious  that  the  depths 
below  the  mixed  layer  are  not  affected  by  the  model. 

The  TOPS  profiles,  on  the  otherhand,  are  affected  at  all 
levels  (down  to  500m).  The  profiles  are  more  complex  and 
harder  to  interpret,  because  of  the  variations  occurring 
over  the  entire  profile.  The  process  of  energy  addition/ 
subtraction  to  the  mixed  layer  is  no  longer  as  obvious  as 
with  a  bulk  model.  Overall,  changes  over  the  13  days, 
however,  clearly  show  a  substantial  increase  in  the  mixed 
layer  temperature;  the  SST  increased  approximately  one 
degree,  while  the  mixed  layer  depth  increased  to  approximately 
40m.  (Figs.  60,  61)  (TOPS  and  TEOTS  profiles  for  the  other 
days  may  be  found  on  the  composite  figures  to  be  discussed 
in  the  last  section.) 

A  comparison  of  the  contours  for  the  TOPS  predicted 
profiles  (Fig.  62)  with  those  of  the  Garwood  model  gives 
further  indications  of  the  difference  in  results  of  the 
diffusion  and  the  ’slab’  or  integral  models.  For  these 
contours  to  resemble  those  obtained  from  the  Garwood  model, 
the  mixed  layers  would  indeed  have  to  be  'slab-like'.  (The 
'plateaus'  located  in  the  center  of  the  contour  field  are  due 
to  repeated  profiles,  i.e.,  profiles  for  6  and  7  August  were 
not  obtained  from  FNOC.) 


Figure  59.  Contoured  Garwood  mode] 
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Figure  60.  TOPS  and  TEOTS  profiles  for  1  August 

1982.  (Temperature  in  degrees  Celsius) 
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e  61.  TOPS  and  TEOTS  profiles  for  13  August 

1982  (Temperature  in  degrees  Celsius) 
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Figure  62.  Contoured  T/Z  field  composed  of  TOPS  profile? 


Both  of  these  sets  of  profiles  and  contours  were  compared 
with  those  obtained  from  daily  averaged  XBT  profiles.  The 
variability  of  the  upper  water  column  as  noted  in  the  latter 
two  sets  of  profiles  was  not  as  apparent  in  the  series  of 
averaged  XBT  profiles,  as  expected.  The  trend  of  the  latter 
XBT ' s, averaged  in  space  and  time,  had  a  tendency  to  wander 

quite  a  bit.  For  example,  the  SST  from  1  August  to  2  August 

jumped  1C  with  a  corresponding  increase  of  0 . 3C  at  400m. 
(Figs.  63,  64)  Similar  changes  occurred  for  other  profiles. 
For  example,  the  13  August  profile  indicated  a  -1.3C  change 
in  SST  and  a  -0.7C  change  in  the  temperature  at  400m  since 
12  August.  (Figs.  65,  66)  (Average  XBT  profiles  for  the 

other  days  may  be  found  on  the  composite  figures  to  be 

discussed  in  the  last  section.)  Mixed  layers  and  their  tem¬ 
poral  variations  are  not  well  represented  by  these  profiles 
due  to  the  spatial  heterogeneity  within  the  study  domain. 

It  takes  a  stretch  of  the  imagination  to  identify  a  mixed 
layer  in  some  of  the  profiles.  Once  identified,  the  mixed 
layers  do  not  tend  to  exhibit  the  characteristics  expected 
for  mixed  layer  evolution,  e.g.,  the  deeper  layers  are  not 
necessarily  cooler.  The  characteristics  seem  to  be  more  of 
a  ’random'  nature,  i.e.  they  are  affected  by  factors  other 
than  those  of  mixing  and  heating. 

The  use  of  average  XBT  profiles  is  somewhat  of  a  contro¬ 
versial  issue.  One  side  holds  that  the  averaging  of  XBT's 
has  some  worth  in  that  it  will  cancel  some  of  the  spatial 
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65.  Average  XBT  and  standard  deviation  plot 
for!2August  1982.  (Temperature  in 
degrees  Celsius) 
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Figure  66.  Average  XBT  and  Standard  deviation 

blot  for  13  August  1982.  (Temperature 
in  degrees  Celsius) 
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and  temporal  effects.  The  other  contends  that  the  averaging 


of  XBT's  destroys  the  vertical  structure,  and,  therefore, 
should  not  be  the  method  by  which  one  obtains  a  single  daily 
profile  with  which  to  compare  the  model  simulations.  The 
latter  believes  that  a  median  XBT  (based  on  various  cri¬ 
teria)  would  be  a  much  better  'measuring  stick'  for  model 
simulation  success.  Having  gone  through  the  abundant  supply 
of  XBT  traces,  comparing  the  averaged  profiles  with  the 
series,  the  author  is  convinced  of  the  latter  viewpoint, 
though  not  contradicting  that  of  the  former.  The  corollary, 
then,  is  that  the  thermal  structure  is  not  robust  over 
modest  distances  of  a  day's  transit  in  a  slow  ship  (ca.  3ms 
It  is  also  suggested  that  greater  care  should  be  taken  in 
grouping  the  XBT  profiles  in  space-time  for  analysis. 
Indicated  on  the  righthand  panels  in  Figs.  67  and  68  are  a 
series  of  hypothetical  XBT  traces  for  two  cases:  one 
containing  XBT's  with  pronounced  jumps  at  the  mixed  layer 
bases,  and  the  other  containing  XBT's  with  gradually  sloping 
mixed  layers.  The  left-hand  panel  of  each  figure  depicts  the 
average  XBT  trace  for  the  series,  as  well  as  a  median  XBT. 
(The  average  profile  is  the  shallower  of  the  two.)  From 
these  it  can  be  seen  that  averaging  distorts  the  vertical 
structure,  making  the  mixed  layer  depth  shallower,  and  the 
temperature  at  the  base  of  the  mixed  layer  colder. 

Median  XBT  profiles  were  almost  as  hard  to  define  as 
were  the  mixed  layer  depths,  especially  over  a  substantial 
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Figure  67.  Averaging  XBT's  with  steep  thermoclines  A)  Right  panel 
XBT’s  to  be  averaged;  B)  Left  Danel:  Average 
(indicated  by  arrow)  and  median  XBT’s 


Figure  68.  Averaging  XBT's  with  less  steep  therm< 
panel:  XBT's  to  be  averaged;  B)  Left 

(indicated  by  arrow)  and  median  XBT’s 


time  and  distance.  XBT  traces  for  a  day  varied  from  having 
'excellent  bulklike '  mixed  layers  to  having  no  mixed  layers. 
The  SST's  did  not  vary  in  any  predictable  pattern,  i.e., 

SST's  of  XBT  traces  without  mixed  layers  were  in  fact  colder 
than  the  SST's  of  the  traces  with  mixed  layers.  The  transit 
of  the  ship  across  numerous  fronts,  jets,  and  eddies  would 
be  one  explanation  for  these  anomalies  and  ambiguities. 

Instead  of  obtaining  a  daily  median  profile,  profiles 
for  the  same  time  each  day  were  used  as  the  'observed' 
reality.  The  time  chosen  was  1800Z  or  1000  local  time  to 
allow  for  the  largest  sample  of  XBT’s,  and  to  avoid  periods 
of  major  diurnal  heating  or  cooling.  Fig.  69  is  the  T/Z 
contour  for  the  XBT ' s  closest  to  1800Z  each  day,  approximately 
within  30  minutes.  There  were  no  XBT's  available  for  the 
designated  time  on  5,  6,  7,  or  8  August.  Therefore,  the 
profiles  for  5,  6  August  were  made  identical  to  that  of  4 
August.  The  profiles  for  7,  8  August  were  made  identical  to 
the  profile  for  9  August.  The  'wavelike'  nature  and  the  two 
events  present  in  the  average  XBT  T/Z  contour  field  at  2,  12 
August  are  also  located  in  the  1800Z  XBT  T/Z  contour  field. 
Associated  with  these  peaks  or  troughs,  one  might  expect  to 
find  a  change  in  the  mixed  layer  due  to  fluctuations  in  the 
whole  water  column.  Such  in  fact  was  the  case  for  mixed 
layers  of  2  and  12  August. 
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Returning  to  the  results  of  the  Garwood  and  TOPS  model 
runs  discussed  thus  far,  one  must  keep  in  mind  a  number  of 
important  points.  Recall  that  the  Garwood  model  simulation 
shown  was  initialized  with  a  profile  from  outside  the  study 
region  and  with  inputs  of  varying  degrees  of  accuracy.  It 
was  not  known  exactly  how  good  the  boundary  conditions  were, 
nor  how  sensitive  that  model  would  be  to  them.  A  sensitivity 
test  was  performed  to  decide  whether  the  Garwood  model  was 
deepening  the  layer  due  to  the  19.5m  FNOC  winds  being  too 
strong.  This  was  done  by  decreasing  the  drag  coefficient 
by  10%,  equivalent  to  decreasing  the  FNOC  winds  to  95%  of 
their  original  magnitude.  As  can  be  seen  from  profiles  for 
1  and  13  August  1982  (Figs.  70,  71),  the  effects  are  not  very 
substantial . 

In  addition,  it  should  be  remembered  that  the  TOPS  model 
was  initialized  and  run  at  a  point  outside  of  the  cruise 
grid.  It  is  not  known  what  values  TOPS  received  as  initial 
or  boundary  conditions,  only  that  it  received  its  atmospheric 
forcing  from  the  atmospheric  primitive  equation  model  and 
NOGAPS,  and  that  eight  XBT’s  were  used  in  the  TEOTS  analysis. 
Finally  a  slightly  warmer  SST  is  expected  for  the  FNOC 
gridpoint,  because  it  is  further  offshore  than  the  study 
domain.  According  to  the  previous  fortnightly  mean, 
however,  the  difference  should  only  be  about  1C. 

To  make  simulations  more  consistent  within  the  domain,  a 
second  run  was  performed  using  XBT  96  to  initialize  the 
Garwood  model.  This  XBT,  deployed  on  4  August,  was  chosen 
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Figure  70.  Daily  average  and  average  tuned  Garwood 
profiles:  1  August  (Temperature  in 

degrees  Celsius) 
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because  its  location  was  resampled  on  10  August  with  the 
launch  of  XBT  216.  (Fig.  72)  (It  was  hoped  that  by  choosing 
XBT  96  to  initialize  the  Garwood  model,  its  mixed  layer 
evolution  could  be  compared  and  evaluated  with  the  data 
obtained  from  XBT  216.  The  initial  SST  used  was  that  of 
XBT  96.  An  initial  mixed  layer  depth  of  20m  was  "assumed".) 

In  addition,  to  avoid  comparing  simulations  of  two 
differently  forced  models,  the  Mellor  level-2.5  model  was 
obtained  and  modified  to  produce  simulations  using  forcing 
identical  to  that  used  in  the  Garwood  model.  Profiles 
resulting  from  the  model  runs  (Figs.  73-85)  along  with  the 
applicable  1800Z  XBT  profiles,  TOPS  predicted  profiles,  ar.d 
daily  averaged  profiles  are  examined.  (Notes:  A)  repeated 
TOPS,  1800Z,  and  average  XBT  profiles  during  in-port  period; 
B)  TOPS  profiles  correspond  to  1600  local  time,  i.e., 
afternoon  heating  exists;  C)  Mellor  profiles  are  not  daily 
averages  . ) 

Based  on  the  overall  change  in  the  profiles,  Garwood  and 
Mellor  profiles  cooled  and  deepened  the  mixed  layer.  The 
profiles  seemed  to  be  quite  consistent  with  the  observed 
data  by  4  August  and  even  more  so  by  10  August.  The  cooling 
trend  in  the  XBT's  seemed  to  be  consistent  with  that  apparent 
from  XBT  96  to  XBT  216.  However,  caution  is  required, 
because  the  role  of  the  horizontal  and  vertical  advection 
has  not  been  assessed. 


150 


CHOSEN  XBT  PROFILES 


TEfffERRTURE 


Figure  72.  XBT’s  96  and  216  (Temperature  in 

degrees  Celsius)  launched  at  approximately 
the  same  location  on  4  and  10  August, 
respectively 
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DAILY  PROFILE,  8/2 
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Figure  74,  August  2  composite  of  Simulated  and 
observed  profiles 


Figure  75.  August  3  composite  of  simulated  and 
observed  profiles 
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DAILY  PROFILE,  8/4 
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Figure  76.  August  4  composite  of  simulated  and 
observed  profiles 
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DAILY  PROFILE,  8/5 
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Figure  77.  August  5  composite  of  simulated  and 
observed  profiles 
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DAILY  PROFILE,  8/6 
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Figure  78.  August  6  composite  of  simulated  and 
observed  profiles 
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DAILY  PROFILE,  8/7 
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Figure  79.  August  7  composite  of  simulated  and 
observed  profiles. 
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DAILY  PROFILE,  8/8 
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August  8  composite  of  simulated  and 
observed  profiles 
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Figure  81.  August  9  composite  of  simulated  and 
observed  profiles 
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DAILY  PROFILE,  8/10 


TEMPERATURE  (C) 


Figure  82.  August  10  composite  of  simulated  and 
observed  profiles 
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DAILY  PROFILE,  8/11 


TEMPERATURE  (C) 


Figure  83.  August  11  composite  of  simulated  and 
Observed  profiles 
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84.  August  12  composite  Of  simulated  and 
Observed  profiles. 
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DAILY  PROFILE,  8/13 
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Figure  85.  August  13  composite  Of  simulated  and 
observed  profiles 
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The  TOPS  profiles,  though  not  expected  to  resolve  the 
small  scale  variability  of  the  boundary  current  region,  on 
the  otherhand,  started  out  warmer  than  the  other  profiles  and 
indicated  additional  warming  with  time.  The  fact  that  the 
gridpoint  used  was  further  offshore  than  the  cruise  region 
caused  one  to  expect  a  slightly  greater  surface  temperature 
for  the  FNOC  profiles.  That  the  forcing  averaged  184  Watts/ 
m**2  per  day  was,  however,  not  anticipated.  This  heating 
is,  however,  believed  to  be  part  of  a  larger  problem,  viz., 
a  heat  bias  in  the  atmospheric  PE  model.  Heats  biases  of 
varying  signs  have  been  previously  identified  by  Elsberry, 
et.al.  (1979),  Budd  (1980),  and  Warn-Varnas,  et.al.  (1982) 
and  shown  to  possess  temporal  and  spatial  scales  of  a  month 
and  tens  of  thousands  of  kilometers,  respectively.  Warn- 
Varnas,  et.al.  (1982),  for  example,  found  that  a  limited- 
area  version  of  TOPS  initialized  and  forced  with  FNOC  fields 
revealed  considerable  skill  during  60-day  runs  in  comparison 
to  persistence  if  a  bias  of  about  150  ly/day  was  removed 
from  the  daily  net  surface  heat  flux  of  the  domain.  This 
correction  proved  to  yield  consistent  results  with  the  heat 
fluxes  computed  directly  using  atmospheric  observations 
acquired  by  ships-of-opportunity .  [Clancy,  et.al.,  1980] 

In  another  study,  comparison  of  -  e  FNOC  heat  fluxes  with 
those  calculated  by  Nate  Clark  indicated  a  cold  bias  believed 
to  be  mainly  the  result  of  excessive  latent  heat  loss.  The 
problem  of  heat  bias  is  extremely  evident  in  such  data  sparse 


areas  as  the  high  latitudes.  For  example,  TEOTS  SST’s  in 
the  Sea  of  Okhotsk,  the  Bering  Sea,  and  the  Labrador  Sea 
have  in  the  past  exhibited  a  warm  bias,  running  as  much  as 
6C  above  those  of  conventional  EOTS  (essentially 
representative  of  climatology  within  these  areas).  [Clancy, 
et.al,  1982] 

In  the  simulations  presented,  the  initial  conditions  and 
forcing  for  the  first  two  days  were  supplied  to  TOPS  by  the 
PE  model.  On  3  August,  NOGAPS  replaced  the  PE  model. 

However,  any  spurious  heat  available  in  the  PE  model  entered 
into  the  simulation  prior  to  bringing  NOGAPS  online.  Once 
introduced  into  the  run,  it  is  unlikely  that  the  heat  could 
be  removed  without  much  mixing  or  surface  cooling  (due  to 
the  Second  Law  of  Thermodynamics).  NOGAPS  *  contribution  to 
the  spurious  heating  is  unknown,  though  it's  clear  that  the 
heating  trend  was  not  corrected  over  the  short  run  time. 

(It  is  as  yet  unknown  whether  NOGAPS  has  a  heat  bias  of  any 
kind. ) 

Examining  the  TEOTS  SST  anomalies  (TEOTS  -  climatology) 
for  some  of  the  cruise  days  as  well  as  for  a  number  of  days 
following  the  cruise,  and  following  the  'return  to  climatology' 
adjustment  installed  in  TOPS  on  20  September  (Table  IV) 
revealed  an  excess  of  2-3C  in  the  study  domain  as  compared 
to  conventional  EOTS. 
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TABLE  IV 


TEMPERATURE  ANOMALIES 


DATE 

TEMPERATURE  ANOMALY  (C) 

820730 

TEOTS  EOTS 

2 

820803 

NOGAPS  Online 

820804 

2  -1 

820808 

2  0 

820811 

3 

820814 

2 

820819 

3 

820827 

2 

820828 

System  Restart 

820910 

1.7 

820920* 

1.4 

820923 

1.9 

♦Return  to  climatology  adjustment  installed 


The  warmer  SST's  along  with  the  an  identified  tendency 
for  the  turbulent  parameter  scheme  used  in  TOPS  to  underpredict 
the  MLD's  [Martin,  1982]  might  be  important  in  explaining  the 
trend  in  the  TOPS'  profiles.  The  difference  between  various 
model  profiles  with  depth  is  believed  to  be  the  result  of  the 
strong  connection  TEOTS  has  with  climatology. 
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To  make  a  strict  comparison  of  model  results  (assuming 
ideally  that  the  model  forcing  is  the  same,  as  it  should  be), 
one  should  interpolate  between  FNOC  gridpoints  to  the 
locations  of  XBT  96. 

It  is  felt  ,  however,  that  the  heating  difference  is  not 
totally  a  function  of  distance.  Erroneous  heat  fluxes, 
strong  dependence  upon  climatology,  and  advection  (associated 
with  mesoscale  variability)  are  thought  to  play  important 
roles  in  the  results. 


V,  SUMMARY  AND  CONCLUSIONS 


The  Garwood  and  Mellor  level-2.5  simulations  resulted  in 
cooling  and  deepenirg  of  the  mixed  layer  over  time,  and 
exhibited  a  likeness  to  the  observations.  The  TOPS  model, 
receiving  different  forcing,  produced  warmer  profiles  with 
a  substantial  warming  trend  over  the  13  days. 

The  exhibited  differences  of  the  TOPS  profiles  with 
respect  to  the  profiles  produced  by  the  other  models  are 
believed  to  be  more  than  a  result  of  location.  They  are 
believed  to  be  due  to  a)  a  heat  bias  in  FNOCrs  atmospheric 
PE  model,  b)  too  much  dependence  upon  climatology,  and 
c)  advection  and  the  subgrid  scale  variability  existing  in 
the  region. 

The  autocorrelation  functions  for  SST  and  MLD  indicated 
.>  a  maximum  correlation  distance  of  35  km.  However,  TOPS 
uses  a  standard  FNOC  grid  with  spacing  381  km  at  60N. 

Hence,  though  interpolation  should  be  carried  out  when  looking 
at  regions  between  the  gridpoints,  it  seems  highly  unlikely 
to  be  able  to  recover  the  small  scale  variability  dominating 
the  region. 
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VI.  RECOMMENDATIONS 


First-  of  all,  numerous  recommendations  are  in  order 
concerning  the  data  acquisition  on  future  cruises.  Most 
importantly,  all  easily  observed  data  should  be  recorded 
regardless  of  the  immediate  mission  of  the  cruise,  (e.g., 
wet  bulb  temperature)  The  equipment  should  be  in  a  working 
status,  and  be  sufficiently  backed-up  with  manual  logs. 
Underway  analysis  should  be  carried  out  in  a  consistent 
manner  wherever  possible  to  avoid  excessive  reworking  of  the 
data  ashore.  A  thorough  scientific  crew  indoctrination 
period  is,  therefore,  in  order. 

From  the  modeling  aspect,  many  more  evaluations  of  TOPS 
are  necessary  under  numerous  varied  conditions  prior  to  its 
full  acceptance  by  the  scientific  community.  It  is  important 
that  the  TOPS  products  be  interpolated  to  the  observed 
position  and  time  of  the  study  domain  prior  to  any  further 
evaluations  of  the  model's  performance.  Design  tests  should 
be  carried  out  with  control  over  the  variables.  Independent, 
direct  measures  of  atmospheric  forcing  should  be  made  to 
drive  the  research  and  operational  models,  (Simulations 
should  be  produced  using  both  sets  of  forcing  values). 
Additionally,  it  would  be  of  interest  to  attempt  to  predict 
SST  and  MLD  from  one  realization  to  the  next,  and  to  examine 
the  changes  in  the  profiles  in  the  (z,x,t)  domain  instead  of 
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merely  the  (z,day)  domain.  All  of  these  tests  and  comparisons 
should  also  be  performed  on  the  60  km  and  20  km  grids  of  the 
future . 

In  addition,  other  operational  models  and  their  products 
are  becoming  more  important  both  in  the  Navy  and  in  the 
civilian  community.  With  their  growing  availability,  so  too 
must  grow  the  amount  and  quality  of  verification  testing 
accorded  them,  for  they  are  an  important  asset  of  the  future. 


1 


J 


171 


APPENDIX  A 


A.  DATA  PROCESSING  PROCEDURES 

Working  with  a  new  data  set  is  rarely  an  easy  task.  The 
data  must  be  edited  and  processed  into  a  usable  format. 

Here,  the  editing  procedures  used  to  produce  the  five 
required  inputs  for  the  version  of  the  Garwood  model  at  hand 
are  described  in  more  detail  to  better  assess  the  model 
inputs  and  results. 

The  procedure  applied  to  the  DAS  data  was  as  follows. 
First,  the  data  on  magnetic  tape  was  dumped  onto  mass  storage 
and  filtered  for  blocks  of  zeros  and  gross  errors  in  lati¬ 
tude  and  longitude.  Following  this,  the  data  from  the 
entrance  and  exit  tracks  to  the  cruise  grid  were  deleted  so 
as  to  confine  the  data  to  the  study  domain.  Finally,  the 
finer  errors  in  temperature,  wind  speed,  etc.  were  deleted. 

The  resulting  air  temperature  data  were  then  averaged  over 
three-hour  intervals  to  yield  the  required  eight  daily  air 
temperatures.  Averages  were  computed  over  bins  containing 
from  6  to  299  values  with  the  maximum  standard  deviation  for 
any  resultant  value  being  1.67C  and  the  minimum  being 
0.0005C.  Only  three  of  the  averaged  values  possessed 
standard  deviations  over  1C.  Procedures  used  to  obtain 
missing  values  within  cruise  days  were  based  on  the 
assumptions  described  previously.  For  days  in  port,  the 
air  temperatures  were  held  at  a  representative  constant  value. 
(This  could  be  improved  upon  with  reasonable  means  to  predict 
gaps . ) 
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The  dewpoint  temperatures  were  erroneous  due  to  equipment 
failure.  The  required  eight  values  of  dewpoint  per  day  were 
calculated  by  assuming  a  relative  humidity  of  90%  and  using 
the  procedure  described  in  Bolton  (1980).  This  value  for 
relative  humidity  was  chosen  to  correspond  with  the 
intermittently  hazy/ foggy  conditions. 

Cloud  cover  indications  from  the  ship's  log  were  either 
in  tenth  of  sky  cover  or  in  descriptive  terms.  The  first 
step  taken  was  to  convert  the  descriptors  to  tenths  of  sky 
cover  using  Bowditch's  American  Practical  Navigator  and  the 
experience  of  the  R/V  ACANIA’s  first  mate.  The  tenths  of 
sky  cover  were  then  converted  to  octals,  and  averaged 
over  three  hour  intervals.  Missing  values  during  cruise 
days  were  not  a  problem  as  they  had  been  with  the  air  and 
dewpoint  temperatures.  The  cloud  cover  during  the  in-port 
days  was  held  at  a  constant  five  tenths. 

Winds  could  have  been  obtained  in  two  ways.  The  DAS 
recorded  winds  could  have  been  converted  to  true  winds  and 
averaged  as  was  done  with  the  air  temperatures.  The  other 
way  winds  could  have  been  obtained  was  by  interpolating  the 
daily  analyzed  TEOTS  39.05N,  127. 75W  wind  to  the  necessary 
three-hour  interval.  The  TEOTS  analyzed  winds  over  both 
CCSII  and  CCSIII  are  depicted  in  Figs.  25  and  26  with  the 
latter  time  series  containing  the  interpolated  eight  three- 
hourly  values.  (The  symbol  denotes  the  value  at  00  Zulu) 

The  latter  method  was  chosen  for  convenience  as  well  as  to 


173 


facilitate  the  comparison  of  the  results  of  the  TOPS,  Mellor 
level-2.5,  and  Garwood  models.  Some  of  the  'CTD  station 
winds’  (already  assumed  to  be  true  winds  due  to  the  'zero' 
speed  of  advance  of  the  ship)  were  plotted  against  the  FNOC 
archived  winds  to  observe  whether  or  not  the  latter  could 
'reasonably'  be  applied  to  the  cruise  area.  The  CTD  winds 
were  found  mainly  to  fluctuate  around  the  FNOC  winds  except 
for  the  first  day  and  a  half  when  the  'CTD  winds’  were 
exceeded  by  the  FNOC  winds.  The  FNOC  winds  were  (Fig.  86) 
employed  in  subsequent  runs.  (The  triangles  represent  'CTD 
winds'  recorded  along  with  suspect  data;  leading  to  a 
possible  suspicion  in  the  winds.)  A  sensitivity  analysis  was 
done  to  assess  the  importance  of  the  wind  stress  as  a  bound¬ 
ary  condition,  and  in  so  doing,  ascertain  the  magnitude  of 
possible  problems  introduced  by  using  winds  recorded  at 
19.5m  vice  10.0m.  The  procedure  and  results  were  described 
in  Chapter  IV. 

The  final  input  required  was  SST.  The  SST's  were  avail¬ 
able  from  a  number  of  sources.  DAS  recorded  boom  temperatures 
could  have  been  averaged.  This  route  was  not  chosen,  because 
of  the  possible  problems  posed  by  the  boom  varying  in  its 
depth  in  the  water  (interaction  with  the  ship's  wake). 
Alternatives  included  the  use  of  the  bucket  or  XBT  surface 
temperatures,  averaged  or  singly,  at  the  required  times. 

For  the  initial  runs,  the  XBT  SST's  recorded  by  DAS 
were  used:  this  was  believed  to  be  the  least  subjective 
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approach.  The  bucket  temperatures  were  subject  to  the 
experience  of  the  watchstander ,  as  were  the  XBT  SST’s 
recorded  in  the  XBT  logs.  The  validity  of  the  latter  were 
found  to  strongly  depend  upon  the  existence  of  a  surface 
thermocline,  either  'real1  or  possibly  due  to  the  initially 
nonlinear  slow  response  rate  of  the  XBT  recorder.  The  XBT 
temperatures  were  taken  from  those  launched  closest  to  the 
required  input  time.  (usually  within  30  minutes — sometimes 
linear  interpolations  were  performed.)  Averages  were  not 
made  for  the  three-hour  transit  of  the  ship  across  possible 
fronts . 

Finally  it  was  necessary  to  input  an  initial  temperature 
profile  into  Garwood  and  Mellor  level-2.5  models.  Temperatures 
were  interpolated  to  intervals  of  one  meter  to  a  depth  of 
200M.  The  first  runs  were  made  using  an  interpolated  TEOTS 
profile  corresponding  to  1  August.  (The  original  TEOTS 
analyzed  temperature  profiles  consisted  of  13  temperature 
depth  values.  Surface  values  were  not  available;  the  tempera¬ 
ture  at  one  meter  was  used  at  the  surface.)  Subsequent  runs 
were  made  using  XBT  96  as  the  initial  temperature  profile  for 
the  Garwood  and  Mellor  level-2.5  models.  This  XBT  was  chosen, 
because  of  its  proximity  to  the  FNOC  gridpoint  of  interest. 

More  importantly,  however,  it  was  chosen  because  of  its  common 
position  with  an  XBT,  216,  launched  six  days  later.  XBT’s  96 
and  216  were  launched  on  4  and  10  August,  respectively,  thus 
spanning  a  part  of  both  cruises. 
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TABLE  V 

VERTICAL  RESOLUTION  OF  TOPS 
VERTICAL  GRID  SPACING  USED  IN  TOPS 


Level 

DepthCm) 

1 

2.5 

2 

7.5 

3 

12.5 

4 

17.5 

5 

25.0 

6 

32.5 

7 

40.0 

8 

50.0 

9 

62.5 

10 

75.0 

11 

100.0 

12 

125.0 

13 

150.0 

14 

200.0 

15 

300.0 

16 

400.0 

17 

500.0 

TOPS  is  run  on  the  Cyber  205  computer  on  a  standard 
63x63  FNOC  global  grid.  Values  of  T,  S,  U,  and  V  are  defined 
with  this  horizontal  resolution  at  each  of  the  seventeen 
levels  listed  above.  All  spatial  derivatives  are  centered 
in  space.  The  horizontal  advection  components  are  defined 
on  a  staggered  grid;  they  are  forward-differenced  in  time. 
Vertical  eddy  fluxes  and  W  are  defined  midway  between  these 
levels.  The  former  are  differenced  backward  in  time,  while 
the  latter  (as  well  as  the  Coriolis  term)  is  differenced 
trapezoidally  in  time. 
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Additionally,  TOPS  is  quasi-three-dimensional  model, 
(though  a  nonadvective  version  is  kept  on  the  Cyber  175  as  a 
backup.)  Advection,  as  previously  mentioned,  is  by 
climatologically  averaged  geostrophic  currents  and 
instantaneous  wind  drift.  The  geostrophic  currents  are 
updated  monthly  using  density  fields  computed  from  monthly 
temperature  and  salinity  climatology  and  the  thermal  wind 
equations  integrated  upward  from  a  predetermined  "level 
of  no  motion”.  The  geostrophic  currents  are  nondivergent . 
The  sole  contributor  to  the  vertical  advection  component 
is,  therefore,  the  instantaneous  wind  drift  via  the  Ekman 
suction  mechanism. 
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